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^ ' Abstract 

I The high- and low-frequency descriptions of the pre-decouphng plasma are deduced from the 

'Tj [ Vlasov-Landau treatment generalized to curved space-times and in the presence of the relativistic 

("^ ' fluctuations of the geometry. It is demonstrated that the interplay between one-fluid and two-fluid 

95 ' treatments is mandatory for a complete and reliable calculation of the polarization observables. 

ff^ ■ The Einstein-Boltzmann hierarchy is generalized to handle the dispersive propagation of the 

■^ ■ electromagnetic disturbances in the pre-decoupling plasma. Given the improved physical and 

^^ . numerical framework, the polarization observables are computed within the magnetized ACDM 

C^ , paradigm (mACDM). In particular, the Faraday-induced B-mode is consistently estimated by 

^ I taking into account the effects of the magnetic fields on the initial conditions of the Boltzmann 

'k^ I hierarchy, on the dynamical equations and on the dispersion relations. The complete calculations 

\-{ ' of the angular power spectra constitutes the first step for the derivation of magnetized maps of 

■ - - ' the CMB temperature and polarization which are here obtained for the first time and within the 

minimal mACDM model. The obtained results set the ground for direct experimental scrutiny of 

large-scale magnetism via the low and high frequency instruments of the Planck explorer satellite. 
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1 Plasma hierarchies and plasma descriptions 

Prior to matter-radiation equality and throughout decouphng the number of charge carriers present 
within the Debye sphere is, overall, inversely proportional to the baryonic concentration. The corre- 
sponding (dimensionless) plasma parameter [H El [3] is O(10~^) and, more specificallja 
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where Xp is the ionization fraction and 
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is the ratio between the baryonic concentration and the photon concentration. Equations (ll.ip - (ll.2p 
assume the values of the cosmological parameters imphed by the 5-year WMAP data alone IH 13 [6|, O [8] 
and analyzed in the light of the conventional ACDM paradigm Q Since the plasma is globally neutral 
(i.e. rie = rii = riQ = ?7bo^7o)) the concentration of charge carriers entering Eq. fll.ip will be ten billion 
times smaller than the photon concentration: this is the ultimate rationale for the hierarchy provided 
by Eq. (II. ip and for the intrinsic validity of the plasma approximation. 

The minuteness of ^fpiasma determines various hierarchies between the physical quantities character- 
izing the pre-decoupling plasma. Every time a hierarchy arises, a potentially interesting (approximate) 
physical description is at our disposal. For instance, the hierarchy between the Hubble rate and the 
collision frequencies of the Thompson and Coulomb scattering permits, before equality, to treat the 
baryon-lepton-photon system as a unique dynamical entity. The latter approximation is always imple- 
mented, in standard Boltzmann solvers, to avoid the stiffness of the numerical system prior to equality. 
In the present paper it will be argued that effective (i.e. one-fluid) descriptions of the baryon-lepton 
system are not adequate for the calculation of the polarization observables when large-scale magnetic 
fields 1^ intervene in the physics of the pre-decoupling plasma. In different words, when the large-scale 
magnetic fields are present it is certainly allowed to solve numerically the system by making use of 
the hierarchy of different physical scales. These hierarchies are, however, less conventional from the 
ones arising in the standard case (i.e. when magnetic fields are absent). The calculation must then 
treat appropriately the specific hierarchies of the problem to avoid that hugely different physical scales 
appear simultaneously in the same numerical integration. In this introduction we first intend to make 
the various hierarchies explicit. In the second place the main strategy of the present calculation will 
be summarized and contrasted with previous attempts. 

To introduce the aspects of the problem in physical terms, it is appropriate to show how the main 
hierarchies of the pre-decoupling plasma are all controlled, directly or indirectly, by fifpiasma- This 



^ Typical plasma parameters in glow discharges are of the order of 10~^. The units used in the present investigation 
will be such that h — c — ks = I- To facilitate the comparison with experimental data, the angular power spectra will 
be however assigned in units of (pK)'^ as it will be explicitly pointed out. 

"^For the explicit estimates, the 5-year WMAP data alone will be consistently adopted. In the numerical study, 
however, different data sets will also be discussed. 

^By large-scales we shall mean here typical length-scales L at least of order of the Hubble radius at equality, i.e. 

rH(Tcq) = H~^\ 



exercise will pin down, implicitly, the different plasma descriptions emerging in the analysis. In short 
the different hierarchies determined by (^piasma are: 

• the hierarchy between the Debye length (i.e. Ad) and the Coulomb mean free path (i.e. Acoui); 

• the hierarchy between the plasma frequency (of the electrons) and the collision frequency; 

• the largeness of the conductivity in units of the plasma frequency. 

To appreciate the validity of this statement, it suffices to write the aforementioned quantities in the 
form of dimensionless ratios depending solely upon Qpiasma- By doing this it can be concluded, for 
instance, that Acoui, is not the shortest scale of the problem: 

Ad fi'plasma , . . 18\/2 

- ^ InAc, Ac = , (1-3) 
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where In Ac is the Coulomb logarithm [21 [3]. With similar manipulations, the plasma frequency of the 
electrons, i.e. Upe, turns out to be much larger than the collision frequency (related to the Coulomb 
rate of interactions), i.e. U 
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Equation (11.41) implies that the electrons will oscillate many times before undergoing a collision. The 
third of the aforementioned hierarchies stipulates that the conductivity (denoted by a) is be paramet- 
rically larger than 1 in units of the plasma frequency: 



a=^E^ = ^^^P^. (1.5) 
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Recalling Eq. (II. ip . Eq. (II. 5p implies, indeed, o jujpc ^ 1 and also that a jT ^ 1. 

According to Eq. (II. ip . the Debye length around equality is of the order of 10 to 100 cm in 
comparison with rn = ii~^ which is, around equality, 20 orders of magnitude larger. Large-scale 
electric (rather than magnetic) fields are highly suppressed by powers of o~^ in the baryon rest frame. 
At the same reference time, magnetic fields can be present over typical length-scales L > L^j where L^. 

T^ / 'r X 1/2 

L. ^ (47rai7eci)-\ ^ = ^T^ — ' (1-6) 
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(see also Eq. (II. 5p ) is called magnetic diffusivity length. For typical values of the cosmological 
parameters, around equahty, L„ ~ 10~^^rH. Magnetic fields over typical length-scales L ~ O(rH) 
(and possibly larger) can be present without suffering appreciable diffusion. Since the magnetic fields 
touched by the present discussion will have, at most, nG strength at the onset of galaxy formation, 
the Larmor radius around equality will be much smaller than the range of variation of the magnetic 
field, i.e. 

TBe < -^ - r-H, r-Be = , V± ~ tJth, (1-7) 



^In this introduction we shall not dwell on the distinction between comoving and physical frequencies. In curved 
space-times of Friedmann-Robertson- Walker type, the electron and proton masses break the conformal invariance and, 
therefore, the plasma and Larmor frequencies will also depend upon the scale factor, as it will be explicitly shown. 



where Vth — JT /rrie and uj-Qe is the Larmor frequency. As pointed out after Eq. fll.6p L > L^ and it 
is always much larger than r-Q^- 

The plasma hierarchies introduced in Eqs. (11. II) . fll.3l) - fll.5p and fll.6p - fll.7l) imply different de- 
scriptions valid for complementary branches of the spectrum of plasma excitations. More specifically: 

• for typical length scales much larger than the Debye length and for typical times exceeding the 
inverse of the plasma frequency, a single-fluid theory naturally emerges and it often dubbed 
magnetohydrodynamics (MHD in what follows); 

• in the opposite limits the two-fluid nature of the system cannot be ignored and, unlike in the 
MHD description, the electromagnetic disturbances propagate in a dispersive medium which has, 
in our case, a finite concentration of charge carriers and a large-scale magnetic field. 

The MHD limit possesses various sub-limits which make the whole dynamics rather rich (see, for 
instance, [H]). The stochastic magnetic field evolving according to MHD and coupled to the fluctuations 
of the geometry affects, by its presence, the propagation of electromagnetic disturbances in the plasma: 
electromagnetic waves with positive (or negative) helicities will propagate with different phase (and 
group) velocities. The latter effect cannot be treated within a one-fluid approximation where the 
displacement current is consistently neglected and the Ohmic current is solenoidal [2]. This is the 
strategy followed in the present paper and it corresponds to the logic used in the analysis of laboratory 
plasmas. The physical and technical challenges of the problem reside in the occurrence that, prior to 
photon decoupling, the metric is a dynamical quantity and that it can fluctuate. It is easily imaginable 
that, in spite of the initial physical analogy with laboratory plasmas, the actual problem will be 
technically more difficult than in Minkowskian space-time. 

Accurate calculations of magnetized CMB anisotropies are an essential tools for the scrutiny of the 
origin and evolution of large-scale magnetism which is observed in the largest gravitationally bound 
systems such as clusters (TUl E], galaxies [121 [12] and even in some superculster. For the interplay 
between CMB physics and large-scale magnetic fields see, for instance, p3l[T5] . More general reviews 
on the problems and challenges of large-scale magnetism can be found in [T6[ [T71 [18] . In simple words, 
we want to embark in these calculations because we want to give reasonable predictions of the potential 
effects of large-scale magnetic fields on the CMB anisotropies. 

Large-scale magnetic fields present prior to equality affect the evolution of the fluctuations of the 
scalar modes of the geometry. The evolution equations of the Boltzmann hierarchy (as well as the 
initial conditions) must be appropriately modifled. This theoretical problem has been scrutinized in 
a number of recent papers [HI [201 [21] (see also [22]). In [19] the estimate of the Sachs- Wolfe plateau 
has been carried on by employing the technique of the transfer matrices. In [20] the temperature 
autocorrelations and the polarization cross-correlations have been computed numerically in the tight- 
coupling approximation. In [21] a semi-analytical evaluation of the TT angular power spectra has been 
carried on[^. 

A dedicated numerical approach for the calculation of magnetized temperature and polarization 
observables has been devised in a series of recent papers [231 IMl |25] . Such a numerical approach in 



^Following the current terminology we will denote by TT the angular power spectrum of the temperature autocor- 
relation. By TE we shall denote the power spectrum of the cross-correlation between temperature and polarization. 
Finally the EE spectrum denotes the autocorrelation of the polarization. The precise definition of the various power 
spectra is discussed in Section [5l 



constructed from one of the standard Boltzmann solvers i.e. CMBFAST [261 l2Z] which is, in turn, 
based on the COSMICS package by E. Bertschinger [281 129] • In [23] the minimal framework for the 
analysis of the effects of large-scale magnetic fields on the CMB anisotropics has been spelled out and 
dubbed magnetized ACDM model (mACDM model in what follows). In the minimal realization of 
the mACDM scenario the inclusion of large-scale magnetic fields amounts to the addition of two extra 
parameters. In [23] it has been shown that the numerical analysis leads to shapes of the TT angular 
power spectra which are exactly the ones computed in [21]. The numerical approach is intrinsically 
more accurate especially at high multipoles. In [25] all the possible (non-adiabatic) initial conditions 
of the magnetized Einstein-Boltzmann hierarchy have been worked out analytically and scrutinized 
numerically. The results reported in [191 ISOl [21] and in [231 1231 [25] bring the treatment of magnetized 
CMB anisotropics to the same standards employed in the case when large-scale magnetic fields are 
absent from the very beginning. 

The problem left out from previous analyses, as stressed in [21], has to do with a more consistent 
calculation of the polarization observables. The problem is, in short, the following. The large-scale 
description of temperature anisotropics demands a coarse grained (one-fluid) approach for the electron- 
ion system: this is the so called baryon fluid which is treated (with no exceptions) as a single fluid in 
popular Boltzmann solvers such as COSMICS, CMBFAST and their descendants. On the other hand 
the dispersive propagation of electromagnetic disturbances demands to treat separately electrons and 
ions, at least at high frequencies. Faraday rotation is one of the situations where the inadequacy of the 
one-fluid approximation (for the baryon- lepton fluid) is manifest. The positive and negative helicities 
composing the (linear) CMB polarization experience, in a background magnetic field, two different 
phase velocities, two different dielectric contants and, ultimately, two different refractive indices. The 
mismatch between the refractive index of the positive and negative helicities induces, effectively, a 
rotation of the CMB polarization and, hence, a B-mode. The inclusion of the Faraday effect in the 
treatment implies, physically, that the proton-electron fluid (sometimes dubbed as baryon fluid) should 
be treated as effectively composed by two different species, i.e. the electrons and the ions. 

Various studies were concerned, in the past, with the estimate of Faraday rotation effects in the 
framework of CMB physics. None of the previous studies, however, could profit of a dedicated nu- 
merical approach where the effects of the magnetic fields could be included at the level of the initial 
conditions and at the level of the dynamical equations. Before outlining the analytical and numerical 
strategies used in the present paper the main results obtained so far will be summarized. 

In [30l [3ll [32] it has been noted, within slightly different perspectives, that the (hnear) CMB 
polarization can be Faraday rotated. A common aspect of the attempts of [501 [SB [22] was that 
the magnetic field was assumed to be uniform (i.e. homogeneous in space) and described within a 
simplified magnetohydrodynamical description which consisted, effectively, to enforce the conservation 
of the magnetic flux. In [32] , on the basis of an explicit model, it was argued that the magnetic fields 
should be treated, in fact, as stochastically distributed if we do not want to break (exphcitly) the 
spatial isotropy of the background geometry: the large-scale magnetic fields arise typically from the 
parametric amplification of vacuum fluctuations or from some phase transition (see, e.g. [T6j) and in 
both cases the produced fields are not uniform but rather stochastically distributed. 

The first convincing measurements of CMB polarization [331 1331 |35] (see also [361 ISZ]) clearly 
suggested the adiabaticity of the initial conditions (because of the location of the first anticorrelation 
peak in the TE power spectra). In [38] the TB correlations have been computed from the initial TE 
correlations provided by the adiabatic mode. The main assumption of [38j has been that the magnetic 



field is, once more, uniform and tliat it does not affect, in any way, tlie initial TE angular power 
spectra. In [39j and [50] it was recognized that the uniform field case is not realistic. The interplay 
between stochastic magnetic fields and Faraday rotation has been more directly scrutinized. The work 
of Ref. [ID] suggested to treat the effect by applying a Faraday screen to the CMB polarization. Also 
in [lO] the main assumption has been to neglect any possible effect of the magnetic fields on the TE 
and EE angular power spectra. 

In the present paper, it will not be assumed that the large-scale magnetic fields are uniform. Hence, 
the direct effects of the magnetic fields on the TE and EE angular power spectra will not be neglected. 
The line of reasoning pursued in the present paper is based on the plasma hierarchies outlined in Eqs. 
(11.11) and in Eqs. f ll.3l) -( TL5l) . The strategy will be to generahze, around matter-radiation equality and 
throughout decoupling, the standard treatment of weakly coupled plasmas in the different branches 
of the spectrum of plasma excitations. The smallness of (^piasma (together with the smallness of the 
electrons and ions kinetic temperatures in comparison with the corresponding masses) allows to enforce 
the (cold) plasma approximation to an excellent degree. The fiat space-equations cannot be simply 
employed around decoupling or even equality since space-time is not fiat. Furthermore, the fluctuations 
of the geometry must be properly taken into account since they are still relativistic at the time when the 
initial conditions of the temperature and polarization anisotropics are customarily set [T9 | [20 t [2T].The 
physical rationale for the strategy employed here is rooted on the simple observation of Eq. (II. 7p . 
which stipulates that the Larmor radius is always smaller than the typical inhomogeneity scale of the 
magnetic field. The approximation which will be adopted here is called, in plasma physics, the guiding 
center approximation and it is due to the pioneering work of Alfven [ITl W2\ . 

A relevant class of results of the present investigation concerns the problem of simulating maps of 
magnetized CMB anisotropics. Magnetized maps of the temperature and polarization observables will 
be reported to illustrate the viability of our numerical approaclo. In the light of the Planck explorer 
mission [43j it will be particularly important to have magnetized maps both for the temperature auto- 
correlations and maps for the polarization observables. This step is mandatory once the corresponding 
angular power spectra can be numerically computed. With our numerical code, appropriately extended 
to include the Faraday mixing term and the two-fluid effects, we are able to compute accurately all 
the required power spectra. Armed with all these necessary theoretical tools, magnetized CMB maps 
can be obtained. 

The present paper is organized as follows. In Section 2 the generalized two-fluid description for the 
charged species (electrons and ions) is introduced. In Section 3 the two-fluid treatment is shown to be 
equivalent to a single fluid description at large-scales. In Sections 4 and 5 the high-frequency branch 
of the spectrum of plasma excitations is discussed. Sections 6 and 7 are devoted to the calculation 
of the TE, EE and BB angular power spectra. Illustrative examples of magnetized CMB maps are 
collected in Section 8. Section 9 contains our concluding remarks. To avoid excessive technicalities, 
relevant results and derivations have been included in the Appendix. 



^Even if temperature autocorrelations are not central to the present discussion (which is focussed on magnetized 
polarization observables) we will also report maps of the temperature anisotropics, mainly for completeness. 



2 Electrons and ions 

There is the custom, in CMB studies, to treat ions and electrons as a single dynamical entity (see, for 
instance, [29] and also [HI HHl |46] ) . This is certainly justified for typical scales L ^ Ad (see also Eq. 
(11.11) ) and for frequencies parametrically smaller than the plasma frequency. Owing to the largeness 
of the Coulomb rate ions and electrons are tightly coupled and the effective degree of freedom which 
should be studied is the so-called baryon velocity, i.e. the centre-of-mass velocity of the electron-ion 
system. Furthermore, owing to the largeness of Thompson scattering, baryons and photons are also 
tightly coupled (but just well before equality). This observation is the basis of the tight-coupling 
expansion which was pioneered in [S] (see also [17]) and since then widely used for semi-analytical 
estimates of the CMB temperature autocorrelations. 

If large-scale magnetic fields are included in the game some of the considerations of the previous 
paragraph apply [19], [20l |2T] (see also [231 121] )■ However, if we ought to treat phenomena related to 
the propagation of electromagnetic disturbances in a plasma, the one-fiuid description is known to be 
insufficient [3]. The insufficiency of the one-fiuid description is already apparent in the calculation of 
the magnetized polarization observables. This caveat has been made explicit in |24 j . 

In what follows the pre-decoupling plasma will be treated as a laboratory plasmeu with the following 
notable differences: 

• around radiation-matter equality and thought photon decoupling space-time is curved; 

• since electrons and ions are non-relativistic conformal invariance is explicitly broken and this 
fact has various physical implications, as we shall see; 

• relativistic fluctuations of the geometry, immaterial for plasmas in the laboratory, have to be 
included in our discussion. 

The conformally flat line element characterizing the space-time around equality and decoupling will 
be written as 

ds^ = g^ydx^dx" = a'^{T)[dT'^ — dx'^], g^y = o?{T)rj^y, (2.1) 

where r is the conformal time coordinate to distinguish it from the cosmic time t (of course dt = 
a{r)dT). The perturbed metric will be parametrized, without loss of generality, in the synchronous 
coordinate system where 13 

SsQijix, t) = a^{T)hij{x, t). (2.2) 

Equation (12.21) holds in real space. It is often useful to separate the fluctuation of the geometry in a 
trace part supplemented by a traceless contribution. For practical reasons this is done in Fourier space 



^It is amusing to notice that the plasma parameter of the decouphng plasma (i.e. 0(10^^), see Eq. (jl.ip ) is of the 
same order of magnitude of the plasma parameter of a tokamak. Of course, in a tokamak, the density of charge carries 
will be much larger. However, the Debye length will be comparatively much smaller. For instance, in a tokamak the 
concentration of charge carriers will be typically 10^'' cm~'^; the typical temperature will be of the order of 10 keV. 
Consequently, the Debye scale is, in this situation, 7 x 10~^ m; the corresponding plasma parameter (defined in Eq. 
(jl.ip ) will then be 0(10^^). This, as anticipated, is the same figure one obtains around equality. 

^° For practical purposes the pivotal description will be taken to be the one of the synchronous gauge [23l EH [25] . 
Different descriptions (in complementary gauges) can be straightforwardly deduced using the fully gauge-invariant 
approaches of [l9l[20] and [87]. 



where the perturbed metric becomes, 



^sgi.jik, t) = a (r) kikjh{k, r) + 6^(A;, r) ( kik 



(2.3) 



which makes clear that the hij{x,T) of Eq. fl2.2p carries, effectively, only two (scalar) degrees of 
freedom. If a fully inhomogeneous magnetic field is present, the evolution equations of the scalar 
inhomogeneities will be necessarily affected. The homogeneous evolution of the background geometry 
(see Eq. (12.11) ) is dictated by the Friedmann-Lemaitre equations: 

37f = 8TTGa^Pu H^ - W = AnGa^pt + pt), p[ + 3n{pt + pt) = 0, (2.4) 

where the prime denotes a derivation with respect to the conformal time coordinate and where Ti = 
a' /a. The quantities pt and pt denote the total energy density and the total pressure of the plasma 
which consists, in the present situation, by two charged species (electrons and ions) and by neutral 
species three neutral species (neutrinos, photons and CDM particles). The dark-energy component 
will be simply parametrized, as in the case of the ACDM scenario, by a cosmological constant. 

The evolution equations of the problem will now be introduced with particular attention to the 
electromagnetic part. There is the custom of treating dispersive phenomena within a single fluid 
plasma description. This habit is, strictly speaking, misleading as already pointed out in ^48j. The 
logic followed in the present treatment can be summarized as follows: 

• start with a bona fide two-fluid description in curved space taking into account all the relevant 
degrees of freedom and, in particular, the electrons and the ions as separated components; 

• derive the one-fluid description which will allow to follow the evolution of the inhomogeneous 
magnetic field; 

• derive the relevant dispersion relations. 

A reasonably general description of the electron-ion system can be achieved from the evolution equation 
of the one-body distribution functions which have the characteristic Vlasov form (see also Appendix 
El for a derivation): 

^ + v-W,f, + eiE + vxB)-V,^f,= (^) , (2.5) 

or \OT J coll 

^ + ^- V.^/e -e{E + vxB)- V,-/e = (1^) ^^, (2.6) 

where the collision terms are provided by Coulomb scattering and where the electric and magnetic 
fields are rescaled as: 

B = a^B, E = a^S, v = ^ (2.7) 

In Eqs. (12. 5p and (12.61) . g is the comoving three-momentum which is customarily defined in the curved- 
space Boltzmann treatment, i.e. q = ap. For a massive particle, such as the electron or the proton, the 
conjugate momenta satisfy Qa/sP'^P^ = Tn?- By identifying g.ijP^P^ = —Sijp'^p^ , the comoving three 
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momentum is introduced I. If we would be in the ultra-relativistic limit (where electrons and protons 
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See Appendix [XI for further details on the geodesies of massive particles endowed with an electric charge. 



are effectively massless) Eqs. (12.51) and (12.61) (written in terms of the conformal time coordinate r and 
in terms of tlie appropriately rescaled fields of Eq. (12.71) ) would be exactly the same equations we 
would have in fiat space-time. This is a consequence of the conformallly fiat nature of the background 
(see Eq. (12.11) ). Prior to equality, i.e. when the calculation of the CMB anisotropies is initialized, the 
electrons and ions are non-relativistic and the presence of a mass term will break (explicitly) conformal 
invariance. This observation is clear by noticing that, in the non-relativistic limit, Eq. (12. 7p implies 
that V = q/{ma). 

It is useful, for immediate convenience, to introduce here the distinction between cold and warm 
pla^mafl The pre-equality plasma is cold in the sense that the kinetic temperatures of the charged 
species are much smaller than the corresponding masses. In the opposite case the plasma is warm. If 
the initial conditions of the Boltzmann hierarchy are set at a temperature which is say one tenth of the 
temperature of neutrino decoupling, the plasma will be already cold to a good approximation. In can 
be shown that, in a warm plasma approach, the corrections to the Faraday rate will be suppressed, to 
leading order, by powers of T^/me [IHl l49] . 

The (approximate) equilibrium solution of the Boltzmann equation will be a Maxwellian velocity 
distribution and, from this observation, the Boltzmann equations can be perturbed to obtain the 
evolution equations of the various moments of the distribution functions such as the evolution of the 
charge concentrations (from the zeroth-order moment), the evolution equations for the velocities (i.e. 
the first-order moment) and so on. While more detailed considerations can be found in Appendix \^ 
the same two-fluid equations obtainable from the Vlasov description can be recovered from the charge 
and four-momentum conservation in curved space-time, as it will be discussed in a moment. 

The evolution equations of the gauge fields will be given, in the present context, by the appropriate 
Maxwell equations which can be written aa^^l 

V.F^'' = 47r(jr + j:), V,F^' = 0, (2.8) 

V^jT = 0, V^j: = 0, (2.9) 

where V^ denotes the covariant derivative. The electron and ion currents are given, respectively, by 

j,^ = -ene<, jr = eni<, (2.10) 

where ni and fze are, respectively, the ion and electron concentrations. The velocity fields satisfy the 
conditions 

g^puy^ = 1, (7,^<Mf = 1, (2.11) 

which imply, in the non-relativistic limit, that Mg = m? = 1/a. The Maxwell field strengths and their 
duals are simply 

Fo, = a^8,, Fi, = -a\,,kB\ F^' = -^, ^^ = ^S^. (2.12) 

Using Eqs. (12.101) . (12. lip and (I2.12p inside Eq. (12. 8p . the following set of equations can be obtained: 

V-^ = 47re(ni-ne), (2.13) 



^^This distinction is customarily employed in plasma literature (see, e. g. [31 [2]) 

^•^The covariant and controvariant indices of the various vectors and tensors must not be confused with the subscripts 
(always in ronian style) which denote the corresponding species. 
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V-B = 0, (2.14) 

VxE + B' = 0, (2.15) 

V X B = AneiriiVi - n^v^) + E' . (2.16) 

where, E = a?£ and B = a?B are the rescaled electric and magnetic fields. The concentrations and 
velocities appearing in the four-currents (see Eq. (12.101) ) have been rescaled as: 

nj = a'^fii, Tie = a? he, 

v\ = aui, Ve = aue- (2-17) 

where v^ and v^ are the conioving three- velocities o In terms of the comoving concentrations defined 
in Eq. (I2.17p . the charge conservation equation, i.e. Eq. (12. 9 p will imply, for electrons and ions, 

dn 



^^ + ^m; + t/i ■ Vn; = 0, (2.18) 



dr 



-F 6'ene + We ■ Vne = 0. (2.19) 



where 6'i = V ■ iTi and ^c = V ■ We are the three-divergences of the comoving three-velocities. Since the 
plasma is globally neutral, the electron and ion concentrations are equal and are approximately ten 
billion times smaller than the photon concentration, i.e. Ue = ni = hq = ?7bo^70 where ?7bo has been 
given by Eq. (11.21) . Owing to this figure the metric fiuctuations can be neglected in Eqs. (12.181) and 
(I2.19P (see, however, Eqs. (JA.ISP and (]A.14p for their explicit inclusion). 

Equations (I2.18P and (I2.19P (derived from Eq. (12. 9p ) can be obtained directly by taking the 
zeroth-order moment of the evolution equations for the one-body distribution functions. In particular, 
defining 



ne(x, t) =no d^vfe{x, V, r), ni(x, t) = hq (fvfi{x, v, r), (2.20) 

Ve{x,T) =no d^Vvfe{x,V,T), Vi{x,T) = Uq I d^v V fi{x,v, t), (2-21) 



we indeed obtain, from Eqs. (12. 5p and (12. 6p . 

_i + V ■ {n,v,) = 0, ^ + ^ ■ Kt/e) = 0, (2.22) 

which are equivalent to Eqs. (l2.18p and (I2.19p . 

The evolution equations for the velocity fields and the density contrasts can also be derived by 
perturbing, to first order and in the synchronous gauge (see Eq. (12. 3p ). the covariant momentum 
conservation: 

9,4T;;f + d^T^^jTi; + r;:„4T^r + sr^pT^^ + r^^^Tjf = f^-j^^\ (2.23) 

9,4Tg) + ^s^U^Z) + K^^^T^'^ + ^^KpTl^ + r^/sTjf = F^'^3^\ (2.24) 

where the barred symbols denote the background values of the corresponding quantity and where ^s 
represents the scalar fiuctuation of a given tensor or connection. Recalling that the electrons and ions 



^"^These variables emerge naturally from the geodesies of charged particles in the gravitational field. They are relevant 
for a correct derivation of the Vlasov-Landau equation in the configuration-velocity space. 



77^00 pe 


77^00 Pi 




r rpOO Pi"i 


S.rpOi _ Pe i 




Pi = "^i^i, 


. <5sPe 



are both non-relativistic, the corresponding energy-momentum tensors will be given as PI 

(2.25) 

(2.26) 
(2.27) 

where 

. <5sPe (5sPi 

Pe = mene, Pi = m^riu o^ = , di = . (2.28) 

Pe Pi 

The masses of the electrons and ions are given, respectively, by rrie = 0.511 MeV and rrii ~ nip = 0.938 
GeV. The binary collision between electrons and protons are rather efficient in bringing the whole 
system to an approximate common temperature which will coincide with the photon temperature 
because of the strength of Thompson scattering (see also p4]). 

Within the same conventions of Eqs. (I2.25p . (I2.26P and (I2.27p . the time- like and space-like compo- 
nents of Eqs. (I2.23P and ( 12.24^ lead, respectively, to the following equations 

S'e = -ee + ^- ^,E . Ve, (2.29) 



h' en-i 



S[ = -0i + - + ^E ■ V;, (2.30) 



CTl -» -* 

V,' + nVo = j[E + Ve X B] + Cep, (2.31) 

PqCl 
Vi' + nv, = ^[E + VixB]+Cpe, (2.32) 

where C^p and Cpc follow from the collision terms provided by Coulomb scattering. Needless to say that 
Eqs. (I2.29p - (l2.30p and 02.31j) - (l2.32j) can also be directly obtained from the moments of the Vlasov 
equation (written in the synchronous gauge) as discussed in Appendix \M 

For typical length-scales much larger than the Debye length and for angular frequencies u -C ujpe, 
Eqs. (^M>-(^M> and flOT]) - fl^:^ reduce to an effective one-fluid theory. In the one-fluid limit, 
by definition, the propagation of the electromagnetic (i.e. high frequency) disturbances is negligible 
and the displacement current vanishes. The Ohmic electric fields are then vanishing in the baryon 
rest frame and the (fully inhomogeneous) magnetic fields affect the Boltzmann hierarchy through the 
magnetic pressure, the magnetic energy density, the anisotropic stress and the Lorentz force. 

In the opposite limit (i.e. angular frequencies larger than the plasma frequency), the two- fluid 
nature of the problem becomes physically relevant: the one-fluid approximation breaks down and 
electromagnetic waves will not simply travel at the speed of light (i.e. 1 in our units). Their group 
velocity will be affected both by the plasma and Larmor frequencies of electrons and ions. Finally, 
since the motion of electrons and ions is non-relativistic (i.e. T -C rrie^p) we also have, in general terms. 



that Xu^pe = \/T/{2me) -C 1. In the following two sections the one-fluid and two- fluid treatments 
will be physically compared in the light of the numerical results which will be later reported. 



^^ Since the electron and ion pressures are given by pc = ficTc and by pi = fiiTi they are suppressed as T/toc,p and 
shall be neglected. 
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3 Large-scale magnetic field and one-fiuid treatment 

Defining, respectively, the centre of mass velocity of the electron and ion system v\y, the total current 
J and the total density contrast (5b 

Vh = ■ ^^, J = e{niVi-neVe), d^ = ■ , (3.1) 

me + mp Pe + Pi 

the generalized MHD reduction can be also implemented in curved space-times. In Eq. ( 13. ip v^ 
is, effectively, the bulk velocity of the plasma: vy, (and its three-divergence 9^) is exactly what it is 
normally employed to describe, in standard Boltzmann solvers [26l [271 EHl EH] , the baryon velocity. 
Conversely (5b denotes the total density contrast of the electron-ion system and it is often naively 
identified with the baryon density contrast. The evolution equation for ^b can be derived by taking 
the three-divergence of the evolution equation of the electrons and of the ions. From Eq. (I2.32p . since 
V ■ -E = we obtain 

P77 — » — * — » 

e', + ne, = — '-v-iv,xB), e, = \/-v,, (3.2) 

Pea* 

e{ + nOi = —^v -{vixB), e, = v- v,. (3.3) 

Summing up Eq. (13.21) (multiplied by rrie) and Eq. (13. 3p (multiplied by nip) the evolution equation for 
^b immediately follows. Similarly, the evolution equation for the total density contrast of the electron 
ion system can be derived. The resulting equations for 6'b and (5b is finally given by: 



9y, + 'Hd^a = 7 s- + - — e {0^ — ^b 



PbaMl + — ^ ^^"^ 



m 



(3.4) 



^b = -^b + ^ + ^^^, (3.5) 



where the Thomposn drag term has been included as: 



e 



e' = neXeCTThO, ^Th = -vtI ^ ) , (3.6) 

where e' is the differential optical depth (xe is the ionization fraction). Note that e' has been correctly 
written in terms of fie (and not in terms of our comoving concentration n^ = a^fie). The Thompson drag 
is dominated by the electron-photon cross-section since the photon-proton cross section is, roughly, 
six orders of magnitude smaller. In Eq. (13. 4p . 6^ is the three- divergence of the photon velocity. The 
photons are neutral species and, therefore, they will be described, in the first approximation, by the 
perturbed covariant conservation equation which leads to: 

0', = -f^% + e'(^b - e,), (5; = -p, + 1. (3.7) 

Well before photon decoupling the photon and baryon velocity fields can be further combined. Indeed, 
during the radiation epoch, the efficiency of the Thompson drag synchronizes the baryon and photon 
velocities almost exactly, i.e. 6*^ ~ 6'b = 6'^b- Indeed, by subtracting Eqs. (13. 4p and (13.70 . the evolution 
equation for {6^ — 6^,) can be obtained and it has the form: 



^/ , /r^b + l\,^ ^ , V-[Jxfi] , V^^^ 



(^b - 0,)' + e' -^ (6^ - e,) = 4 i + ^ - ^^b, (3.8) 
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where the rrie/mp <^ 1 has been consistently neglected and where 

3pb /690.18\/ uj^o \ 

is the baryon to photon ratio. The ratio m^/m-p is customarily used as an expansion parameter in 
the standard MHD approximation scheme |2]. Equation (13.81) shows that, in spite of the pre-equality 
differences in 9^ and ^b the two velocities will be inevitably driven towards a common value, i.e. ^-yb 
which follows from the sum of Eqs. (13.41) and (13. 7p : 

where p = a^p^. While the sum of the momentum conservation for the ions and for the electrons led 
to the bulk velocity of the plasma, another linear combination of the same set of equations leads to 
the Ohm's law, i.e. 



e^no 



rcoui/= 0, (3.ii: 



E + v{,y. B 
.rrieaj L 

where the Hall terms (proportional to J x B) has been neglected and where the collision terms have 
been estimated in the linear approximation around a Maxwellian distribution of velocities. The neglect 
of the Hall term is perfectly consistent since we will keep only terms which are quadratic in the magnetic 
field intensities. This term would lead in the Einstein equations to terms that are cubic in the magnetic 
field intensities and shall then be negligible. Equation (13. lip leads then to the canonical form of the 
Ohmic current o 

J = a{E + v^xB), a=^^ = -.f^. (3.12) 

It can be easily checked that the expression of the conductivity of Eq. (I3.12p is exactly the one 
anticipated in Eq. (II. 5p . In this one-fluid limit the displacement current vanishes and, therefore, 

J=— Vxfi, E + v^xB = — , 3.13 

47r 47ro" 

which demonstrates that large-scale electric flelds are suppressed in the baryon rest frame by one 

power of the conductivity. Furthermore, collisions are sufficiently frequent to keep the isotropy and 

this fit with the stochastic nature of the pre-decoupling magnetic fleld. Equation (I3.13P can be used 

to simplify the evolution equations of the photon-baryon velocity. In particular, using simple vector 

identities it can be easily shown that the following relations are satisfled by the large-scale magnetic 

fields: 

J-V-[/xfi] = VVB-^, (3.14) 

where ^b and a-Q are two dimensionless quantities which are related, respectively, to the magnetic 
energy density and to the anisotropic stress. Because of the large value of the conductivity, the electric 
component of the MHD energy-momentum tensor are suppressed by two powers of the conductivity: 

5X = 5sPB. 47;^ = -4pB + fl|, 5,T^ = ^\^^f/ (3.15) 

4Pj3 = _^, 4P3 = ^, Ul = J_(BiB^--6iBA, (3.16) 



^^For some applications is useful to recall that a — aa where a is the physical conductivity and a is the comoving 
conductivity. The expression of a is obtained from a by writing T = aT, i.e. by going back from the comoving to the 
physical temperatures. 
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where B^ = BiB^. Defining then 

Qb = ^, djd'fll = {P, + P7) Wb, (3.17) 

P7 

Equation (13 .14^ can be recovered by repeated use of known vector identities. Note that the Poynting 
vector of Eq. (I3.15P is only suppressed by one power of the conductivity since it contains the Ohmic 
electric field which vanishes as cr~^ in the baryon rest frame. 

The MHD description allows then to reduce the problem of the evolution of large-scale magnetic 
fields to a rather well defined system where the plasma is globally neutral and the Ohmic current is 
solenoidal. This description is ideal in order to study the large-scale effects of the magnetic fields, i.e. 
exactly when gravity is important. Indeed the MHD energy- momentum tensor must be consistently 
included in the perturbed Einstein equations which take then the form: 

2 V'e + 'Hh' = -SirGa^ [S^Pt + ^sPb] , (3.18) 

V^^' = 4na\pt + pt)eu (3.19) 

h" + 2nh' + 2V^e = 247rGa^[6sPt + Spb], (3.20) 

{h + 60" + 2n{h + 60' + 2V'e = 2ATTGa^[{p, + p,)a, + (p^ + p^)(Tb], (3.21) 

where Eqs. (I3.18P and (I3.19P are, respectively, the Hamiltonian and the momentum constraints. 
Equations (I3.20p and (I3.2ip follow, respectively, from the {i = j) and (i 7^ j) components of the 
perturbed Einstein equations. The term containing au is simply related to the neutrino anisotropic 
stress as did^U) = {j)u + Pv)V^(Ty. In Eqs. (I3.18p - (l3.2ip d^Pt and 6't denote the total density fluctuation 



and the total velocity field: 

4pt = 4P7 + ^sPu + 4pi + (5pc + (5pc, (3.22) 

4 4 

(Pt + Pt)^t = -P7^7 + -pyO^ + Pe6'e + pA + PcOc- (3.23) 

In Eqs. (I3.22P and (I3.23P the various subscripts refer, with obvious notation, to the various species of 
the plasma and the different numerical weights in Eq. (13.230 are a simple reflection of the difference of 
the various pressures. The only species not introduced so far are the CDM particles and the neutrinos. 
The CDM component obeys, in the synchronous gauge, the following pair of equations 

<5: = -^c + y, e'^ + m, = Q. (3.24) 

The neutrino component decouples at temperatures of the order of the MeV and it is therefore de- 
scribed by the appropriate Boltzmann equation which reads, in the synchronous gauge (see last part 
of Appendix |X]) , 

Tl + ikpTy = -40 + 2p^{h' + 60). (3.25) 

The evolution equations of lowest moments of the Boltzmann equation can be obtained with standard 
techniques and they are: 

5'. = -p. + Ih', (3.26) 

C = '^Su - A;V„ (3.27) 

4 SA- 2 

^:. = Y^^.-^-^.3-Y^(/^' + 60), (3.28) 
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where a^ = Tyil'^-, i-e. the neutrino anisotropic stress appearing in Eq. (13.211) is simply related by a 
numerical factor to the quadrupole moment of the neutrino phase space distribution. 

Most of the considerations related to the presence of the magnetic fields before equality can be 
conducted in rather general terms. For instance, the analytical form of the initial conditions (reported 
for completeness in Appendix [B]) can be deduced without specifying the magnetic field configuration. 
Prior to equality magnetic fields can certainly be present even if they are physically constrained by 
the isotropy of the background geometry: the magnetic field cannot be uniform. If the magnetic field 
is stochastically distributed, its two-point function in Fourier space must be divergenceless and it is 
given by 

kB,i$)B,m = ^-^nAk)VB{k)6^'\k + p), P,,{k) = (6., - ^), (3.29) 

where Ps(/c) denotes the magnetic power spectrum. In principle there could be a second contribution 
to the two-point function of Eq. (I3.29p . Such a contribution would lead to a non- vanishing magnetic 
helicity. Configurations of this type might arise at the electroweak epoch and have been named 
hypermagnetic knots [SHI EI] (see also [521 ES])- As far as the scalar modes are concerned the hehcal 
part of the field plays no role. The same conclusion holds in the case of the Faraday effect by a 
stochastic magnetic field, as it has been explicitly discussed in [39] (see also [lO]). 

It is useful to mention, at this point, that large-scale (tangled) magnetic fields might have also 
specific effects related to the vector and tensor modes of the geometry (which are minute at large 
scales). These effects have been analyzed in [3HIS51ISB]- The analysis of large-scale magnetic fields can 
also be conducted within fully covariant approaches [STJ [58] which are, when appropriately handled, 
equivalent to the one devised in this paper and based on earlier analyses [191 



4 High-frequency waves and two-fluid effects 

As already noticed in [2l| , the tight Thompson coupling, well verified prior to equality, effectively breaks 
down, by definition, around photon decoupling so the baryon and photon velocities (i.e. 9^ 7^ ^7) will 
have to be numerically followed, as in [21|, by means of a generalized slip equation. In similar terms, 
when the frequency of the polarized CMB photons greatly exceeds the plasma frequency, we can expect 
that the two-fluid nature of the baryon-lepton system will be somehow "resolved" . By resolved we 
mean that it will make actually a difference that instead of one (globally neutral) fluid the plasma is 
formed by two (intrinsically charged) fluids (i.e. 9i ^ 9^ and Ve 7^ Vi). 

Dispersive phenomena in the pre-decoupling plasma arise exactly in this limit where the angular 
frequencies of the propagating electromagnetic waves are larger than the plasma frequency. Faraday 
rotation (as well as other dispersive phenomena at high frequencies) can be derived by studying 
the evolution of electromagnetic disturbances when the electromagnetic background contains a finite 
concentration of charges and a stochastic magnetic field. 

This analysis mirrors what is customarily done in the case of a magnetized plasmas in the labora- 
tory: MHD is used to describe the large-scale magnetic field while the two fluid treatment is enforced 
to scrutinize the propagation of electromagnetic disturbances. In this sense, the large-scale magnetic 
field plays effectively the role of a background field. This terminology is a bit ambiguous in the sense 
that the magnetic field is not part of the geometric background (i.e. it does not determine the expansion 
rate at equality or decoupling), but it is rather part of the electromagnetic background determining 
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the dispersion relations for the waves propagating in the plasma. 

Equations (l2A8l) and (l2A9ll together with Equations (1X321) and Eqs. (l233ll - ( 12161 ) can be linearized 
in the presence of the magnetic field B{x), i.e. 

rae,i(T, x) = riQ + 5rae,i(T, x), B{t, x) = B{x) + b{T,x), 

We,i(T, x) = 6vc^i{T,x), E{t,x) = e(r,f). (4.1) 

where B{x) is the large-scale magnetic obeying the MHD equations discussed in the previous section 
and characterized by the two-point function given in Eq. ( 13.291) . The field variables e and b denote, 
respectively, the electric and magnetic fields of the wave. The absence of a background electric field is a 
reflection of the largeness of the Coulomb conductivity in units of the plasma frequency (see Eq. (11.51) 
and also Section [3]). Using the notations of Eq. (14. ip . Eqs. (12.181) . (12.191) and (I2.32p lead, respectively, 
to the following set of equations: 

(5< + rioV ■ 6ve = 0, 6n[ + uqW ■ 5v; = 0, (4.2) 

e 



Sv'^ + nsve 



ruea I 



e + 5ve X B 



5vl + 7i5vi 



e + SifiX B 



(4.3) 



rnia 

Vxe = -6', V-e = 0, (4.4) 

V X 6 = e' + 47reno(5i7i -5i4). (4.5) 

As the structure of the system shows, the waves parallel (i.e. ||) and orthogonal to the magnetic field B 
(i.e. _L)will then obey different dispersion relations. The important proviso is that the magnetic field 
is inhomogeneous over typical scales that are much larger than the Larmor radius. This observation, 
already mentioned in Eq. (11.71) . will now be made more explicit. Let us consider, for instance, the 
decoupling time which takes place, according to the WMAP 5-year data, for Zdec = 1089.9. At this 
epoch Tdcc — 0.25 eV and rn^Tdcc) = 286 Mpc. The typical scale of variation of the magnetic field at 
decoupling, i.e. L(rdec) will be of the order of the Hubble radius and certainly larger than the magnetic 
diffusivity length defined in Eq. ( 11.61) : 



rBe(rdcc) = f^ ^ 1.19 X 10^*^ —— cm, 

L^Tdec) = (47r(Tdeci^dcc)"'/' ^ 1-4 X lo" cm. (4.6) 

These figures imply that, indeed, L(rdec) > ^Be(Tdec) since L(rdcc) ^ -^^(Tdec)- The magnetic field, 
inhomogeneous over scales much larger than the Larmor radius, can therefore be considered uniform 
for the purpose of deriving the dispersion relations. What has been introduced here in terms of the 
typical scales of our problem is indeed one of the most useful approximations of the physics of cold 
plasmas and it has been pioneered by Alfven ^T[ B2] . 

In general the background magnetic fields will be, as in our case, solution of the Maxwell's equations 
and will therefore be both space and time dependent. Let us restrict ourselves, for the present, 
to spatially inhomogeneous magnetic fielddlZl B{x). If the inhomogeneity is small, it is possible to 
determine the trajectory of the charged species as a perturbation around the center of the particle 
orbit: 

B{x) = B{xo) + {6x ■ V)B\s=^-;^, 6x = x- fo, (4.7) 



^^This is the physical situation for B{x) if the magnetic flux is strictly conserved. In other words, it can be shown 
that B{x) is a solution of the MHD equations in the limit of infinite conductivity. In the opposite case a power series 
in l/cr will arise in the background solution. 
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where Xq is the instantaneous position of the guiding center. Equation (14. 7p holds provided 6B (i.e. 
the change in the magnetic field over a typical distance of the order of the Larnior radius r^e) is such 
that \SB\ = \{Sx ■ V)-B| <C B. But this implies exactly Eq. (14. 6p . i.e. ree <^ L where L is the distance 
over which the magnetic field changes significantly. As already mentioned this perturbative approach 
was extensively applied by Alfven and it often dubbed as the guiding centre approximation in the 
plasma literature. This point is rather important and it is precisely where our treatment diverges 
from the studies reported in Ref . [SOt 131] [32| [38] . In those treatments the magnetic fields have been 
considered just globally uniform. This would break spatial isotropy and it is not true for a stochastic 
magnetic field. The correct physical statement is that the magnetic field is uniform over the typical 
scale of the particle orbit, i.e. the field experienced by the electron in traversing a Larmor orbit is 
almost constant. It is then practical to Fourier transform with respect to the conformal time variable 


b{x, t) = b^(x)e-'I^''-, e{x, r) = e^(x)e-^I^^-, (4.8) 

The solution of Eq. (14. 3 p is trivial along the direction of B and the total current can be written as 

where a = iTi/oJ = H/oj <^ 1 accounts for the curved-space corrections and where 



UJr 



AnnQXce"^ _ 
\ rUpa \ 



AirniXee'^ 



nip 



(4.10) 



_ AnrioXce'^ AnneXee"^ 

V rnptt y rup 



In Eqs. (14.1 op and (14. lip the comoving (angular) frequencies uJpe and uJbs have been related to the 
corresponding physical frequencies Upe and c^Be by recalling that no = a'^Ue = a^n-i. It is relevant to 
mention, once more, that the presence of the scale factor in cJpi and tJpe is a direct consequence of the 
breaking of conformal invariance induced by the masses of the electrons and ions. In the direction 
orthogonal to B, the evolution equations of the electromagnetic waves can be recast in the following 
handy form: 

(V X b)zj = -iiJe{uJ, a)ezj, (V x e)zj - iUkj = 0, (4.12) 

where the dielectric tensor e{uJ, a) can be written in a generalized matrix notation as: 

/ e±i(uJ,a) 2e_L2(t^,a) \ 

e(p,a) = -ie±2{uJ,a) e±i(ZJ,a) , (4-13) 

V e\\{uJ,a) J 

having defined the corresponding components as 

e\0, a) = l- J/' , - J^""" , , (4.14) 

e±2(^, a) = ^=2? — r-{\^ — =2 —=2? — r~{\2 — =2~- (^-^S) 



^^In the present and in the following sections w and k denote, respectively, the comoving angular frequency and the 
comoving wavenumber. 
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In Eqs. (14.151) and (14.161) . on top of the comoving plasma frequencies introduced in Eqs. (I4.10p - (l4.1ip . 
there appear also the comoving Larmor frequencies for the electrons and for the ions: 

_ eB ■ h eB ■ h 

UJBe = = ^BeO, Ubs = (4- 17) 

TUea rrie 

_ eB ■ h eB ■ n 

LJBi = = t^BiO, u;Bi = • (4.18) 

rriia Trip 

where the relation between the comoving and the physical magnetic fields, i.e. B = a^B, has been 
exploited. Equations (I4.12p and (14.130 can be finally combined leading to the compact relation 

k'%j, - (k ■ e^,^^)k = eiuJ, a)uJ^e^,-. (4.19) 

Defining 

k\\ = k(k ■ h) = kcos^, k± = ksm{}, (4.20) 

the dispersion relations can be read-off from Eq. (14.190 . in matrix fornj^l: 




[1 - W -^^ 

i^ [cos^t?-^] - sin ?9 cost? | €^,^,^2 | = 0. (4.21) 

— sin -(9 cos -(9 [sin^f? — -^J 

Requiring that the determinant of the matrix (I4.2ip vanishes we get what is sometimes called Appleton- 
Hartree dispersion relation [3]: 

2e|| cos'^e[{n'^ - e-){n^ - e+)] = sin^ 0(e|| - n'^)[n\e+ + e_) - 2e+e_], (4.22) 

where 

e+{uj,a) = e±i{u,a) + e±2{^,(y), e-{u,a) = e±i{u,a) - ei_2{uj,a). (4.23) 

If ^ = in Eq. (I4.22p (i.e. the wave propagates along the magnetic field) the waves with positive helicity 
(i.e. e+) and negative helicity (i.e. e_) experience two different phase velocities given, respectively, by 

v±{uj,a) = — — — -, n±{uJ,a) = Je±{u,a). (4.24) 

n±[uj, ex) 

If the propagation of the wave is orthogonal to the magnetic field direction, there are two possible 
modes the ordinary and the extraordinary wave with dispersion relations given, respectively, by 

2/- N /- N 2r- ^ 2e+(p,a)e-{uj,a) 

n {uj,a] = e\\(u,a], n(uj,a) = — —^ — r —^ — r-. (4.25) 

e+[uj,a) + e_[uj,a) 

The dispersion relations of Eq. ( J4.25P are simply aJ^ = A;^ in the physical range of parameters. The 
reason for this conclusion is that, in the physical system under study, 

• the plasma and ion frequencies of the electrons are always much larger than the corresponding 
frequencies of the ions; 

• the CMB angular frequency is, in turn, much larger than, both, the Larmor and plasma frequen- 
cies of the electrons. 



^^ As usual the refractive index is defined as uj/k = 1/n. The refractive index n should not be confused, of course, 
with the unit vector n. 

17 



The comoving Larmor and plasma frequencies for electrons and ions are, respectively: 

a;Be = 0.01759(— — I Hz, cJei = 9.578 x lO^^f —— IHz, (4.26) 

V nG / V nG / 

aJpe = 0.285 v/^ (3^) '^' MHz, ojpi = 6.652 J^, (3^) ^'^ kHz. (4.27) 

P ^ V0.02773/ ' pi V e ^Q_Q2273; ^ ' 

Since in our code the present value of the scale factor is normalized to 1, the present value of the 
physical frequency coincides, by definition, with the comoving frequencies. The figures of Eqs. (14.261) 
and (I4.27P should be compared with the typical frequency of CMB photons which is clearly much larger: 
the maximum of the CMB emission is located for a comoving angular frequency Umax = 27rz7max where 
i^max = 222.617 GHz. We are therefore in the physical regime specified by the following hierarchies 
between frequencies 

?^»1, ?HL»i, §^»i, ^»i, (4.28) 

i^Be l^Bi ^pi i-t^pe 

where aJ denotes a typical CMB angular frequency. Under the conditions expressed by Eq. (14.281) the 
expressions of e± (cJ, a) greatly simplifies and the result is, quite straightforwardly, 

6+(U,«) = 1 - TlU- V ^^-29) 

a;[a;(a + 1) + c^BeJ 

e_(cJ,a) = l-^- ^^^ .. (4.30) 

Equations (14.29^ and (I4.30p imply, together with the hierarchies deduced in Eq. (I4.28p . that the 
dispersion relations of the ordinary and extraordinary wave are u"^ = k"^. Indeed, the relation between 
the comoving wavenumber and the comoving angular frequency for the ordinary (i.e. ko) and for the 
extraordinary (i.e. /ce) waves): 



Co = ^A 



^pe 



\ uJ\l + a) 



u 



u'(a + iy-uJU-2uLuJ\a + l)+u;l 



''\ uj'[uj\a + ir-uJl^-uJlia + l)] 



(4.31) 



(4.32) 



It is clear that Eq. (I4.3ip leads immediately to fco = ^ since, according to Eqs. (I4.27P and (I4.28p . 
(Upg/cJ^) ~ O{10^^^). Similar conclusion can be drawn, after some algebra, for the extraordinary wave. 
Equations (I4.30p and (I4.32p do not seem to forbid the existence of resonances. The resonance of Eq. 
(I4.30p is the well known cyclotron resonance which is here avoided since the typical CMB frequencies 
are much larger than ubb- For the same reason also the resonance arising in Eq. (I4.32p is avoided for 
the physical range of frequencies. 



5 CMB polarization and Faraday screening 

The different dispersive behaviour of the positive and negative helicities implies a rotation of the 
components of the electric (or magnetic) field of the wave. Thus, also the related Stokes parameters 
will be rotated. Within the adopted set of conventions, the Faraday-rotated Stokes parameters are 
the same ones we would obtain by the appropriate rotation of the coordinate system. Consider, for 



sake of concreteness, a monochromatic wave propagating along the magnetic field direction an linearly 
polarized along ei at r = and z = 0: 



e{z, t) = Eoiie 

Since the positive and negative helicities are defined as 

ei + ie2 



i[u)T—kz) 



ei - 162 



(5.i: 



(5.2) 



72 ' v^ 

the linear polarization is simply composed of two circularly polarized waves, one with positive helicity 
(propagating with wavenumber /c+ = ^je+iuj^ a) oJ) the other with negative helicity (propagating with 

wavenumber k^ = 
helicities: 



e-{uj,a) uj). Equation (15. ip can be rephrased in terms of the two propagating 



eiz^T 



V2 



g g-«(^-r-fc+^) _|_ g_Q-H^T-k+z) 



2 



ei e 



Ak4 



+ e'^'']+ie2[e'^+' 



Ak- 



'^ (5.3) 



where the second equality follows from the definitions of e±; the refractive indices n-|-(cJ, a) = \/t±[uj, a) 
have been introduced in Eq. (I4.24p (see also Eqs. (I4.29p and (I4.30p ). Since the polarization plane of 
the incoming wave is rotated, two out of four Stokes parameter 



/ = |e ■ eip + |e ■ 62^, 
Q = |e ■ eip - |e ■ 62]^, 
t/ = 2Re[(e-ei)*(e-e2)], 
\/ = 2Im[(e-ei)*(e-e2)]. 



(5.4) 
(5.5) 
(5.6) 
(5.7) 



will be rotated: while / and V will be left invariant, Q and U are rotated. Suppose then the initial 
wave is linearly polarized along ei and Q^™-* = Eq. Inserting Eq. (15.30 into Eqs. (15. 5p and (15.60 . the 
Faraday rotated Stokes parameters are 



)(m) 



Qi^) = Qi'-)cos(2A(^^ 



f/(F) = -QM sin (2A^^ 



where 



^'/'^^^ = 77 \/e+(^,a) - ^Je-{u,a) 



Az. 



(5. 



(5.9) 



Equation (15.80 implies also that Q^^yU^^^ = — cot (2Aip'^^^). In more general terms, if the initial wave 
is not polarized along a specific Cartesian direction, Eq. (15.80 becomes: 



Q(F) = gM cos (2A(/?(^)) + f/('°) sin (2A(/?(^)), 

f/(F) _ _g(m) gi^ (2A^(^)) + f/('°) COS (2A<^(^)) 



(5.10) 



In similar terms, while I and V are invariant under rotations of the polarization plane, Q and U do 
rotate if ei and 62 rotate. In fact, defining (^ = Alp^^\ the rotation of the two unit vectors ei and 62 
by ip 

ei ' = ei cos(y9 + 62 siny), 62' = — ei sin;/) + 62 cosy?, (5-11) 

leads to the same rotated Stokes parameters which have been computed in Eq. (I5.10p . 

The rate of rotation per unit time is called Faraday rotation rate and it is given by: 



^(n) 



dr 



2" 



19 



(5.12) 



where we used that, in our units, dz = dr. To compute the difference of the two refraction indices 
appearing in Eq. fl5.12p we can expand e-|-(cZ;, a) given in Eqs. (14.291) and fl4.30p for |cJ/cZ;pe| -C 1 and 
for lu/aJBel < 1- As already discussed (see Eqs. (14.261) . (14.271) and (14.281) ) this is fully justified by the 
physical values of the aforementioned angular frequencies. The Faraday rotation rate then becomes: 



UJBc (UJ 



jr^f,^ = Z^IZ^ 



LU 



B-n 



2Tim1 



~ H/l tpnAJ p 



V 



(5.13) 



where we used that no = a^'n^ and that uj = 27rz/. Equation (I5.13P can be simplified even further by 
taking into account Eq. (13.60 and the definition of the differential optical depth, i.e. e' = aneXeCTh'- 



J^{n) = e'F(n), F(n) 



3 n-B 



(5.14) 



The numerical code used for the calculation is the extended version of the code used in [23l [2ll [25] 
and it will be dubbed, in the following two sections, as MAGcmb. As already acknowledged in the 
introduction, MAGcmb is based on CMBFAST [26l [27] (which is, in turn, based on COSMICS). We 
shall not dwell here on the problem of the initial conditions of the Einstein-Boltzmann hierarchy which 
has been throughly discussed in Refs. [IHl [20] (see also [22]) and in [2ll [25]. By initial conditions 
of the Einstein-Boltzmann hierarchy we simply mean a consistent solution of the Einstein equations 
(see Eqs. (I3.18p - (l3.2ip ) and of the lowest multipoles of the evolution equations of neutrinos, photons, 
baryons and CDM particles in the approximation of tight Coulomb and Thompson scattering. Initial 
conditions are set well before matter-radiation equality. To avoid possible confusions and to make the 
present script self-consistent, the analytic form of the magnetized adiabatic mode has been reported 
in Appendix [B] (see Eqs. (lB.ip - (]B.10p ). Of course, other initial conditions are, in principle, at 
our disposal: they include the magnetized version of all the non-adiabatic modes which have been 
introduced in [20l [211 l22] and discussed, within different perspectives, in [871 EH]. For purposes of 
illustration we will stick here on the minimal mACDM model where the initial conditions are given in 
terms of a single magnetized adiabatic mode. 

The polarization is only generated very near the surface of last scattering as the photons begin to 
decouple from the electrons and generate a quadrupole moment through free-streaming. The effects 
of the magnetic field on the generation of the linear polarization will be first computed according to 
the following set of equations derived in [20l [2l] : 



A;^ + tkfiAQ = e' 
A'u + z/c/iAu = - 



e-^ih' + 6C) 



-Ai + Aio + /iWb - -P2(/i)5'Q 



-Aq + -(1-P2(/^))5q 



e'A, 



< + Hv,, + ^(3zAii + v^) + ik^-^ = 0, 

-Kb 4_Kb 



(5.15) 

(5.16) 
(5.17) 
(5.18) 



where i?b has been defined previously and where we defined Wb = 6\^/{ik). Moreover, in Eqs. (15.150 
and ([516]) : 

5q = Ai2 + Aqo + Aq2. (5.19) 



The notations A^ and Aq£ denote the £-th multipole of Ai and Aq. In Eqs. (I5.15P and (15.161) 
-P2(/^) = (3yU^ — l)/2 is the second Legendre polynomial. 
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From Eqs. fIS.lOp and (15.121) we can take the total time derivative supposing that the initial 
polarization is independent on time. From the Faraday rotation rate we then get: 

A'q + n%Aci = 2e'F(n) Au, (5.20) 

A'u + n'diAu = -2e'F(n)AQ. (5.21) 

Equations (I5.20p and (I5.2ip describe the Faraday rotation mixing. Previous approaches to Faraday 
rotation, consisted in neglecting the effects of the magnetic fields both on the initial conditions of the 
Einstein-Boltzmann hierarchy and on the the dynamical equations. The calculation of the Faraday 
rotation assumed then, as initial condition, just the conventional adiabatic initial condition of the 
ACDM mode. The EE polarization correlations had then no trace of the magnetized contribution. 
It will be shown that the present numerical results, as already suggested in [24J, do not support the 
latter conclusion. 

Before plunging into the presentation of the numerical results obtained within the improved version 
of MAGcmb, two comments are in order. The first comment is merely technical and it has to do with the 
fact that, as it is well known, the transport equations including Faraday terms as well as other dispersive 
effects have been extensively studied in second half of the past century (see, for instance, [521 EDI EI] 
and [62] and references therein). These studies were motivated by highly relativistic astrophysical 
plasmas leading to synchrotron emission which could be influenced, under some circumstances, by the 
Faraday effect. Consequently, the magnetoactive plasma of [59l [60] is relativistic. Furthermore, the 
transport equations do not include the effect of the metric inhomogeneities and are always defined in 
flat space-times, as in the case of the conventional (polarized) heat transfer equations. Finally, with 
few exceptions, the magnetic field is always taken to be uniform or, at most, non- uniform along a 
specific Cartesian direction (see second paper in [60]). 

The second issue we want to discuss has to do with the recent 5-year data release of the WMAP 
collaboration [7] . The WMAP 5- year data have been analyzed to look for possible birefringent effects 
in the polarization observables. We will here scrutinize the parametrization of birefringent effects 
employed in [7] and we will contrast it with the birefringent effects typical of Faraday rotation by a 
stochastic magnetic field. This discussion will also be useful to set up precisely the relation of the 
fluctuations of brightness perturbations in terms of the well known E-modes and B-modes. Two linear 
combinations of the brightness perturbations can be usefully introduced, namely: 

M±{n)=AQ{h)±iA^{n). (5.22) 

Under a rotation of the coordinate system such as the one introduced in Eq. (15. lip we clearly have 
that 

M±{n) = e^^'^M±{n), (5.23) 

where the tilded quantity denotes the transformed combination. Since /iA± transforms as a spin-2 
field, we can expand it in spin-2 spherical harmonics [631 [M], i-6. 



M±{h) = Y^ a±2, em ±2Yem{n), (5.24) 



im 



where ±2yem{n) are the spin-2 spherical harmonics which can be introduced, formally, as Wigner 
matrix elements [HS]- A typical quantum mechanical problem is to look for the representations of the 
operator specifying three-dimensional rotations, i.e. R; this problem is usually approached within the 
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so-called Wigner matrix elements, i.e. T>^^,{R) = {j, m'\R\j, m) where j denotes the eigenvalue of 
J^ and m denotes the eigenvalue of J^. Now, if we replace m' -^ — s, j — > i, we have the definition of 
spin-s spherical harmonics in terms of the 'D_l„i{a, P, 0), i.e. 



syim{a,P) = y^^2)i'U(a,A0), (5.25) 

where a, (3 and 7 (set to zero in the above definition) are the Euler angles defined as in [66]. If s = 0, 
^o,m('^5/^) 0) = \/(2£ + l)/4vrY£m(«,/?) where Yirn{ot-,l3) are the ordinary spherical harmonics. As 
discussed in [65] the spin-s spherical harmonics can be obtained from the spin-0 spherical harmonics 
thanks to the action of certain ladder operators. These ladder operators can be connected to the 
ladder operators of a putative 0(4) group [^. The reason for the appearance of 0(4) stems from 
the fact that we are here composing two rotations: the rotations in the three-dimensional space and 
the rotations around a given point on the tangent plane of the celestial sphere. A rather productive 
approach to obtain explicit expressions for the spin-2 spherical harmonics consists in studying directly 
the polarization on the 2-sphere [HZ]. Either with the ladder operator of [65j or with the more intrinsic 
approach of [67] it is possible to obtain more explicit expressions of the spin-2 spherical harmonics in 
terms of the spin-0 spherical harmonics and the result is: 



±2 



Y,m{^,(l)) 



(^-2)! 



sm V 



y,^(^?,0). (5.26) 



The "electric" and "magnetic" components of polarization are eigenstates of parity and may be defined, 
from a±, nm oi Eq. (15.241) . as 

(E) _ 1/ , N (B) _ "^ / \ (r nn'\ 

^Im — ~-^\0'2,im-r 0,-2,£m), (^tm — 7:\(^2/m. — (l-1,imj ■ l^-^'J 

Under parity inversion, the components appearing in Eqs. (I5.27P and 

aL - (-1)' «L, «L - (-ir^ «F^- (5.28) 

Therefore, the 15-modes have the same parity of the temperature correlations which have, in turn, 
the same parity of conventional spherical harmonics, i.e. (—1)^. On the contrary, the B-modes have 
(—1)^+^ parity. The existence of linear polarization allows for 6 different cross power spectra to be 
determined, in principle, from data that measure the full temperature and polarization anisotropy 
information. The cross power spectra can be defined in terms of the spectral functions C} where 
X and Y stand for E, B or T depending on the cross-correlation one is interested in: 

^|,.,, tj^. (6.29) 



fjiXY) 



Therefore, if we are interested in the TT correlations we just have to set X = T and F = T and use 
the relevant expansions given above. In the present paper, following a consolidated convention, the 
correlations will be denoted, with shorthand notation, as TT, EE, BB, TB and so on. Suppose we are 
interested in the TT correlations, i.e. the usual and well known temperature correlations. From Eq. 
(I5.29P we will have 

c-r' -i-Jk'ikj: fcwi;<"^) . (5.30) 
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Now, using the orthogonality of spherical harmonics, we have that 

<^Jmik) = [ dnY;jn)A,{k,n). (S-Sl^ 



Inserting Eq. ( ]5.3ip into Eq. (]5.30p and recalling the Rayleigh expansion of the fluctuations of the 
intensity: 

Ai(fc, n, r) = E(-0'(2^ + 1) Ak(^, r)Peifi) (5.32) 



we obtain 



^(TT) _ ^ 
71 



fdkP\Au\'^. (5.33) 



To get to Eq. (I5.33P the following two identities have been used, i.e. 

An 
dhPi'ik ■ n)Y;„Xh) = ——-Y;^{k)6ee', 

Y: YLik)Y,Un) = -^Pe{k ■ n). (5.34) 

It will now be shown how the angular power spectrum of the B-mode can be related, with semi- 
analytical methods, to the spectrum of the Faraday rotation rate. At the same time it is rather relevant 
to mention that the birefringence induced by stochastic magnetic fields differs from a particular kind 
of birefringence which has been considered in [7] . When a pseudoscalar field is coupled to F^^^F^^ 
birefringent effects are expected (see, for instance, [IH])- Furthermore, the presence of pseudo-scalar 
interactions was invoked in [SOj [51] (see also [52l [53]) as a mechanism for the generation of magnetic 
helicity from a stochastic background of hypercharge fields with vanishing helicity. Suppose, therefore, 
a coupling in the Lagrangian density of the form (SaF^'^Fap/M where M is a typical mass scale, /5 is a 
dimensionless coupling and a is the dynamical pseudo-scalar field. Let us suppose also, for simplicity, 
that the background magnetic field is absent. This is, in fact, the situation contemplated in [7]. In 
this set-up the evolution equation for the magnetic field of the wave is given by: 

h" - V%+ — a'V X 6 = 0, (5.35) 

Recalling that the prime denotes a derivation with respect to the conformal time coordinate, the shift 
in the polarization plane will be Ax = /5Acr/(2 M) where we denoted by x the angle to avoid confusion 
with the other angular variables previously introduced. Now, it is clear that while Ax is frequency 
independent, Aip^^'^ depends upon the frequency. Furthermore, the inhomogeneity scale of the magnetic 
field plays enters directly Alp^^\ On the contrary Ax = /5Acr/(2M) is a fully homogeneous quantity. 
Now, in [7] Ax is observationally constrained by adding one single parameter to the ACDM model. 

The rotation by Ax implies, at the level of a\^ and of a\J^ 

afl = cos (2Ax)aS + sin (2Ax)aS, (5.36) 

~a2 = - sin (2Ax)aS + cos (2Ax)aS. (5.37) 

where the tilded quantity denotes, as in Eq. 05.231) the rotated variable. Using now Eq. (15.291) 
and assuming that Ax does not depend upon the Fourier wavenumber, the transformed EE and BB 
autocorrelations can be swiftly obtained 

Cr^ = cos^ (2Ax)Cr^ + sin^ (2Ax)Cr\ (5.38) 

Cf ^) = cos' {2Ax)Cf''^ + sin' (2Ax)Cf ""^ (5.39) 

Cr^ = liCr^ - Cr^) sm (4Ax). (5.40) 
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According to the analysis of [7] the tilded quantities denote the observed power spectra. In deriving 
Eqs. (15.381) . (15.391) and (15.401) , the initial C) is assumed to be vanishing for consistency with ^. 
Furthermore, if we would take literally the ACDM model with no tensors, then C) =0 (as also 
assumed in [7]). In the latter case the induced BB angular power spectrum will be simply given by 
Eq. (Km with Cf ^) = 0. 

In the case of birefringence induced by a stochastic magnetic field the calculation is more involved 
and it is reported, in full detail, in Appendix O Using the shorthand notation F{n) = AB ■ h, F{n) 
can be expandes in series of (scalar) spherical harmonics. Thus, since the magnetic field does not 
break spatial isotropy, the angular power spectrum of Faraday rotation can be written as: 

{F{n,)F{n2)) = ^ ^(2^ + ^)CPPi{n, ■ ns) (5.41) 

where 

CP = 47rA'i{i + 1) / Y'^Bik)^^^. (5.42) 

As it will be shown in a moment, the magnetic field affects directly the EE angular power spectra 
and, therefore, the equation for the BB angular power spectrum which follows from Eq. (1C.13I) can 
be written, formally, as 

Cr= E^(^'^i'^2)Cjf^cf (5.43) 

where Q{i,ii,i2) is a rather cumbersome function of the multipole moments which contains also a 
Clebsch- Gordon coefficient. The expression for Q{i, ^1,^2) is reported in Appendix [Uj and has been 
derived in BOI. 



Even if the cases described by Eqs. (I5.38l) -( l5.40l) and Eqs. (I5.42p - (l5.43p are conceptually very 
different, it is tempting to interpret the bounds on Ax as constraints on the averaged rotation angle, 

(iA.<^'R^i:^c.r-/^cf'f. (5.44) 

In spite of the fact that this identification is rather gross and, as a consequence, the derived bounds are 
not so significant and inherently ambiguous. One of the sources of ambiguity stems from the frequency 
dependence of C) which scales, according to our conventions, as (^'max/'^)^ where z/max = 222.617 GHz 
is the maximum of the CMB blackbody o 



The WMAP experiment observes the microwave sky in five frequency channels ranging from 
23 GHz to 94 GHz. The bandwidth, correspondingly, increases from small to high frequencies signal- 
ing that probably the best sensitivity to polarization comes the high frequency channela^^l. On the 
contrary, the Faraday rotation signal is larger at low frequencies and, therefore, a possibility would 
be to compare the bounds on Ax derived in |7j with Faraday rotation signal evaluated at the lowest 
frequency of observation. The lowest frequency of observation [37J is, however, used as a foreground 
template. Therefore the lowest available effective frequency would be indeed around 30 GHz which is, 
approximately, the lowest frequency channel of the Planck experiment [l3j (see also later, in the first 
part of section [7j). The results reported in [7] do not allow to determine how the various frequency 



^°It is practical to refer the frequency dependence to t'max defined as in Eq. (|C.20p and (|C.2ip . Consequently, the 
power spectrum of Faraday rotation will scale as (i^max/^^)'* (see also Eqs. (|C.24p and (|C.25p of Appendix \X^ . 

^^More precisely the five frequency channels of the WMAP experiment are centered at 23, 33, 41, 61 and 94 in units 
of GHz. 
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channels have been combined in the analysis. The authors of [^ simply assumed that the rotation 
angle was frequency independent, which is not our case, as it will also be discussed later on. 



6 Magnetized TE and EE angular power spectra 

Three fiducial sets of parameters will be used for the illustration of the numerical results. The first 
set of parameters is obtained from the WMAP 5-year data alone [H |5l El [7l [8] analyzed in the light of 
the conventional ACDM modern: 



(fibo, ^co, ^A, hQ, ris, r) = (0.0441, 0.214, 0.742, 0.719, 0.963, 0.087). (6.1) 

If compared to the best fit parameters of the WMAP 3-year data alone [36l [37], we can notice: a slight 
increase of the scalar spectral index (from 0.958 to 0.963); a slight increase of fibo; a slight increase of 
r^co and a consequent decrease of ^a. The features discussed with one set of data have a corresponding 
qualitative counterpart within different data sets. The WMAP 5-year data have been also combined 
with the ACBARo data set [68]. In the latter case, the derived values of the cosmological parameters 
are: 

(l^bo, ^co, ^A, /io, «s, r) = (0.0441, 0.215, 0.741, 0.720, 0.964, 0.088). (6.2) 

Finally, a third set of cosmological parameters will be the one suggested by the analysis of the WMAP 
5- year data in combination with the large-scale structure data (in particular 2df and SDDS data) [69j, 
in turn combined with all the available supernova datacj [TOl [71] : 

(fibo, fico, ^A, ho, n„ t) = (0.0462, 0.233, 0.721, 0.701, 0.960, 0.084). (6.3) 

In Fig. [1] the results of the numerical integration of Eqs. (I5.15I) - (I5.18I) are reported for the choice of 
parameters given in Eq. (16.11) . The magnetic spectral index has been fixed, in both plots, to ne = 1.5. 
The purpose of Fig. [T]is twofold. In the first place we point out, as already discussed in ^24j, that we 
are sensitive to nG magnetic fields. The dot-dashed curve in Fig. [1] indeed corresponds to a magnetic 
fielco of 10 nG while the dashed curve corresponds to a magnetic field of 5 nG. 

The second occurrence which is evident from Fig [1] is that the large-scale magnetic fields do affect 
the TE and EE angular power spectra. This cannot happen, by construction, if the magnetic field 
is uniform as in the case discussed, for instance, in [30|, [32] and [38]. Furthermore, for a stochastic 
magnetic field the TB and EB correlations will be vanishing. The same situation illustrated in Fig. [T] 
is also illustrated in Figs. [2] and [3] but with two different sets of parameters. In Fig. [2] the WMAP 
data are combined with the ACBAR data and the full line is then the best-fit with parameters fixed 
as in Eq. (16. 2p . The range of the mACDM parameters illustrated in Figs. [2] and [3] is narrower than 



^^In the present script the differential optical depth has been denoted by e' while a more standard notation is r'. Here 
the variable r denotes instead the conformal time coordinate. To avoid possible confusions in the identifications of the 
mACDM parameters the conventional notation will be here restored. 

^•^The Arcminute Cosmology Bolometer Array Receiver (ACBAR) operates at three frequencies, i.e. 150, 219 and 274 
GHz. 

^''in the plots we will denote this class of data by WMAP 5-year + BAO + SNall. The acronym BAO refers to the 
baryon acoustic oscillations which could be deduced from the large-scale structure data. 

^^The regularized amplitude of the magnetic field, Bl, is defined as in [24]. For blue magnetic spectral index (i.e. 
riB > 1) the window function is taken to be Gaussian. For red spectral index (i.e. jt-b < 1) the window function is 
a step with smooth edges. The magnetic pivot scale, denoted by k^ is always taken, for illustrative purposes, to be 
lMpc~^. In Appendix [Cl the precise definition of Bj^ is reviewed since it is needed for a consistent evaluation of the 
Faraday rotation power spectrum. 
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WMAP 5-year alone, n = 0.963, n = 1 .5 
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Figure 1: The TE and EE angular power spectra in the mACDM model. The full line corresponds to 
the best fit of the WMAP 5- year alone data analyzed in the light of the conventional ACDM paradigm. 

WMAP 5-year -i- ACBAR, n^ = 0.964, B,^ = 5 nG 





mo 


WMAP 5- 


-year+ACBAR, 


n 

s 


=0.964, B^= 5 nG 






A 








' 








5-year best fit 
WMAP +ACBAR 








/ \ 








VO-6 




M 


50 



•. / 


A 


/ 


..^ 


- - V-8 




5~ 

1 


'\ 


1 A 


/ 




-50 


\/ 


M 


/ 


\j y 




-100 




V 








Figure 2: As in Fig. [T] the TE and EE angular power spectra for a different choice of magnetic field 
parameters. With the full line the best fit to the conventional ACDM paradigm is illustrated and it 
corresponds to the parameters reported in Eq. (16.21) . 



in the case of Fig. [TJ Indeed, as it can be argued by comparing the legends of the three figures, 
in Fig. [T] a magnetic field of 10 nG has been allowed. On the contrary in Figs. [2] and [3] a sizably 
smaller magnetic field has been assumed (i.e. -Bl = 5 nG) but the magnetic spectral index ne has 
been illustrated for two different blue values. Different sets of parameters have been studied within 
the mACDM model and the qualitative conclusions are the same. Contrary to what normally assumed 
large-scale magnetic fields do affect both the TE and EE angular power spectra and nG magnetic fields 
lead to clearly observable features. 

To conclude the present section it is relevant to point out, very swiftly, that the TT angular power 
spectra have been studied, both, analytically and numerically in [201 EI] and in [231 121] ■ The typical 
patterns of correlated distortions of the acoustic peaks have been thoroughly investigated (in the light 
of the minimal mACDM model) in [231 121] • In Fig. H] the plots for the TT angular power spectra 
are reported. In the plot at the left the diamonds are the WMAP 3-year data the boxes the WMAP 
5-year data. With the dashed line the best fit to the WMAP 3-year data alone is illustrated. With 
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Figure 3: The same as in Fig. [2] but with the ACDM parameters reported in Eq. ([3]). 
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Figure 4: The magnetized TT angular power spectra in hnear scale (plot at the left) and in double 
logarithmic scale (plot at the right). 

the full line the magnetic filed is included with the other parameters fixed to the 5-year best fit. In 
the plot at the right the TT and EE angular power spectra are reported in a double logarithmic scale 
to show the relative magnitude of the two angular power spectra. 



7 Faraday screened BB angular power spectra 

The C) depend not only upon the parameters of the magnetized background but also upon the 
comoving frequency. This important occurrence, as explained in the previous sections, is a direct 
consequence of the dispersion relations. The Planck experiment [13] will observe the microwave sky 
in nine frequency channels: three frequency channels (i.e. u = 30, 44, 70 GHz) belong to the low 
frequency instrument (LFI); six channels (i.e. u = 100, 143, 217, 353, 545, 857 GHz) belong to the 
high frequency instrument (HFI). The BB power spectra for all the relevant frequency channels. There 
are reasons to expect that the sensitivity to polarization will be larger at high frequency [13]. At the 
same time the expected signal will be larger at small frequencies. Since our source of information is 
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Figure 5: The BB angular power spectrum in linear scale (plot at the left) and in logarithmic scale 
(plot at the right). 



primarily the Planck blueboolo we will illustrate our results for a putative frequency of 100 GHz. 
Needless to say that all the relevant magnetized polarization observables can be computed, with 
MAGcmb, for any relevant frequency and explicit examples will be given in what follows. In Fig. \5\ 
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Figure 6: The effect of the increase in the spectral index in the region of the peaks. 

the BB angular power spectrum is reported for the parameters listed in the title of the figure. The 
difference between the left and the right plot is just the scale which is linear (at the left) and which is 
logarithmic (at the right). While the linear scale is rather good in capturing the region of the peaks, 
the logarithmic scale accounts better for the low multipoles. In Fig. [5l with the full line, the BB 
angular power spectrum is reported in the case when the the EE and the TE spectra are computed 
taking into account the magnetized contribution. The dashed line, on the contrary, is obtained when 
the EE power spectrum is just taken to be the one implied by the 5-year WMAP data analyzed in 
the light of the conventional ACDM model. Figure [5] demonstrates that if the E-mode is computed 
consistently, the resulting B-mode is not the Faraday screened E-mode produced by the conventional 
adiabatic mode. 



^^We are aware of larger sensitivities possibly achievable in the context of the LFI. However, since these results have 
not been published, we will stick to the informations contained in the bluebook. 
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By increasing the magnetic spectral index n-Q the peaks become less pronounced while the corre- 
sponding amplitude of the BB angular power spectrum decreases. This aspect is captured by Fig. [61 
At the left the region of the peaks is illustrated for i?L = 5 nG and ns = 1.5. At the right, for the 
same value of the regularized amplitude, the spectral index is increased from 1.5 to 2. This trend is 
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Figure 7: BB power spectra for different values of the parameters of the magnetized background (plot 
at the left). EE and BB power spectra are compared (plot at the right) with a double logarithmic 
scale. 

even sharper in Fig. [71 In the plot at the left, with the full line, the case of -Bl = 1 nG is illustrated. 
The other three curves (i.e. dot-dashed, dashed and dotted) all refer to the case of -Bl = 5nG but with 
progressively increasing spectral indices. Always in Fig. [71 (plot at the right) the relative amplitude of 
the EE and BB power spectra is illustrated. 

By looking at the relative amplitude of the EE and BB angular power spectra (see, for instance, 
right plot in Fig. [7]), it could be argued that, in principle, the full line (i.e. the BB angular power 
spectrum) and the dashed line (i.e. the EE angular power spectrum) might intersect for frequencies 
^ < i^max- Indeed, while the EE spectrum does not depend upon the frequency, the BB spectrum 
increases as the frequency decreases well below z/max- The analytical and numerical considerations 
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Figure 8: Bounds stemming from the consistency of the calculation, 
pursued so far hold under the assumption of sufficiently small Faraday rotation rate, i.e. F{n) < 1. 
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This is not a severe limitation and does not help to exclude the mutual cross-over of the EE and BB 
power spectra. Indeed, the requirement F{n) < 1 implies a condition on Eq. (15.441) . In particular, 
to keep the estimate semi-analytical, we can observe that the integrand appearing at the right hand 
side always increases (as a function of £) if we consider, just for illustrative purposes, the case ub > 1. 
Thus the integral will be dominated, in the first approximation, by £ ~ imax and this observation 
leads to a (not so constraining) bound which is reported in Fig. [8] (plot at the left). In the {Bi^,n-Q) 
plane the constraints will be different depending upon the specific frequency. Always in Fig. [8] (see 
the plot at the right) the mean square rotation angle (i.e. 



|A<^| = ^(|A^(^)|2), see also Eq. I^M)) is 
reported for different values of the regularized magnetic field intensity i?L and for different frequencies. 
The maximal frequency will be taken to correspond to the frequency of the CMB maximum (i.e. 
i^max = 222.617). The minimal frequency will be taken to be z/min = 30GHz. Now, the WMAP 5- year 
data give a series of bounds on the rotation angle due to the birefringent nature of the primeval plasma. 
As already discussed the analysis of [i7j does not apply to the present case for, at least, two reasons. 
In the case studied in the present paper the rotation angle has a well defined power spectrum (i.e. 
an explicit i dependence). Secondly, the rotation angle also depend upon the frequency. In Fig. M 
with the full thick line we report the bound |A^| < 1, as implied by the internal consistency of the 
calculation. With the dashed (thick) line the bound \A^\ < 0.04 is reported. The latter figure would 
read, in degrees, \A^\ < 2.4 deg and it is one of the most constraining bounds reported in [7]. 

The WMAP 5-year data (see, in particular, [8]) imply that, when averaged over i = 2 — 6, i{i + 
l)Cf, /{2tc) < 0.15(;uK)2 ( 95% C.L.). Slightly less stringent bounds were obtained by using the 
3- year data [37J. As already mentioned the putative constraint of [8J does not make reference to a 
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Figure 9: Direct constraints on the BB angular power spectrum are illustrated, for different frequencies, 
with the full, dashed and dot-dashed horizontal (thick) lines. 



specific frequency. So we have to assume that it holds for all the frequencies. Then, in our case it 
should be imposed at the lowest frequency which is the most constraining one. The lowest available 
frequency for this purpose would be for 27 GHz. The preceding frequency (i.e. 23 GHz) has been used 
as a foreground template [37] and, consequently, the EE and BB power spectra have not been freed 
from the foreground contamination. We therefore choose to set the bound for a minimal frequency of 
30 GHz since this not only intermediate between the KKa and KQ bands of the WMAP experiment 
but it is also the putative (lowest) frequency of the Planck experiment. As expected, the direct bound 
on the BB angular power spectrum is not so constraining. This aspect is visible in Fig. [9] where, for 
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Figure 10: Direct constraints on the distortion and increase of the first acoustic peak of the TT angular 
power spectrum are reported. 

illustration, the case of blue spectrum has been considered. The different bounds (horizontal lines) 
refer to different fiducial frequencies. The full, dashed and dot dashed lines refer, respectively, to the 
cases of z/ = 30 GHz, u = 100 GHz and z/ = i^max- 

In [23l [21] it was argued that an idea on the range of variation of the parameters of the magnetized 
background can be obtained by monitoring the height and the position of the acoustic peaks. In 
particular, the WMAP collaboration measures that the Doppler multipole is given by [7| 

m + i) 



2n 



-a 



(TT) 
Doppler 



5756 ± 42 {fiKY 



(7.1) 



In Fig. [To] the Doppler multipole is computed for different values of the parameters of the magnetized 
background. The full (thick) lines give the region allowed by error bar quoted in Eq. (17.11) . 

It is finally useful to recall that in [2l| it has been suggested that the highest peaks in the acoustic 
oscillations might be sensitive to large-scale magnetic fields. Since the 5-year WMAP data have been 
combined with ACBAR data we might suspect that the bounds stemming from the acoustic peaks 
of higher order (i.e. compatible with the ACBAR range) would put more severe constraints. This 
is certainly a perspective which should be investigated. From our preliminary investigation it seems 
that the height of the first acoustic peak is still one of the best indicators of the typical patterns 
of correlated distortions induced by a large-scale magnetic field. This conclusion can be drawn by 
comparing, for instance, the values of the parameters in Fig. [10] with the plot at the left in Fig. [11] 
A magnetic field of 5 nG and relatively high spectral index is barely tolerable by the ACBAR data 
but it is definitely ruled out by the accuracy on the first peak. Already with a magnetic field of 1 
nG the situation much less constrained. This also shows that our calculation is sensitive to magnetic 
fields of nG strength and smaller. This was never achieved in the past. Needless to say that the high-£ 
region is expected to be accurately measured by the Planck explorer satellite [13] and this was also 
one of the motivations of our analysis. We want to stress that it would be very interesting to achieve a 
reasonable accuracy also on the EE power spectra at high multipoles. In Fig. [11] (plot at the right) we 
report Z)^ = |C^ — C^ \/Ci , where C^ is computed without any magnetized contribution 
and for the best -fit to the 5-year WMAP data supplemented by the ACBAR data. On the contrary 
C) is computed in the presence of a magnetized background with parameters specified in the title 
of the figure. This gives an idea of the required accuracy of the observations at high multipoles. 
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Figure 11: Angular power spectra at high 



Up to now it has been shown how the magnetized polarization observables can be accurately 
computed. This analysis concludes the first part of the program commenced in [23] and further 
developed in [211 125] • The results obtained in the present paper allow to compute not only the TT 
but also the TE, the EE and the BB power spectra which are the ones expected in the case of a fully 
inhomogeneous magnetic fielccj. The next step is to feed MAGcmb in one of the standard strategies 
of parameter estimation and this analysis is already in progress. 

In the case of birefringent effects due to pseudo-scalar interactions the strategy adopted in [7] has 
been to check for a potential reparametrization of the BB and EE angular power spectra. Such a 
reparametrization would be directly induced by birefringent effects and has been also discussed in Eqs. 
fl5.38p and (15.391) . The advantage of such a strategy is that, of course, only one extra parameter is 
introduced on top of the parameter content of the conventional ACDM paradigm. 

In what follows we would like to suggest a complementary strategy which would allow, in principle, 
to determine if the birefringent effects are really due to pseudo-scalar particles or rather to a stochas- 
tically distributed large-scale magnetic fields. The key observation, in this respect, is the following. In 
spite of the fact that the magnetic field intensity enters, in a rather simple way, the Faraday rotation 
rate, the BB angular power spectrum does not have a simple scaling property as a function of -Bl the 
reason is that, as explicitly demonstrated, the magnetic field modifies also the EE and TE spectra 
which are needed, in turn, to compute the BB power spectrum. Conversely, a simple scaling law of 
the signal can be deduced by looking at the frequency dependence. In short the logic is that while 
the EE and the TT autocorrelations do not depend upon the frequency, the BB autocorrelations do 
depend upon z/ as a consequence of the dispersion relations. Denoting by 



,(T) 



+ i: 



27r 



c, 



(TT) 



'(E) 



+ i: 



27r 



■C. 



(EE) 



'(B) 



27r 



-C, 



(BB) 



(7.2) 



the results obtained in the present paper suggest that when the frequency is rescaled a.s u ^> u, the 
three relevant autocorrelations of Eq. (17.21) change as: 
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(7.3) 



^^As already stressed in the introduction, the uniform field case, breaking explicitly spatial isotropy, can only be 
considered as useful toy model which is however unrealistic to begin with. 
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There are two possible objections to this argument. The first one is that the frequency channels should 
be sufficiently numerous. The second objection is instead that the scahng law expressed by Eq. (17.31) 
could be contaminated by spurious scalings induced by the know (and yet unknown) foregrounds. 

It seems that nine frequency channels (such as the ones of the Planck satellite) are sufficient for 
testing the scaling law of Eq. (17. 3p . There is, however, a proviso. It seems that for an effective 
determination of the scaling law (17.31) the polarization sensitivity of the three low frequency channels 
is an essential requirement. 

Usually, in CMB studies, the frequency dependence is not used, as we suggest, to deduce the 
nature of the signal but rather in order to disentangle the possible foreground contaminations. Indeed 
both the synchrotron and the free-free emissions lead to specific frequency slopes cj. So the question 
becomes: are the known foregrounds able to simulate the scaling law of Eq. (17. 3p . The answer to the 
latter question is negative. 

To demonstrate the previous statement let us adopt the parametrization of the foregrounds adopted 
in [37], namely, in the notations of Eq. (17. 2p 
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(7.4) 
(7.5) 



Equation (\7A\\ is a simplified parametrization of the EE foreground while Eq. (17.51) is a parametrization 
of the BB foreground. Both parametrization hold outside the galactic mask (conventionally called 
P06). In Eqs. (17.41) and (17.51) the two components of the foreground are given, respectively, by the 
dust (characterized by the subscript d) and by the synchrotron (characterized by the subscript s). The 
values of the various parameters appearing in Eqs. (17. 4p and (17. 5p are: 

QsE = 0.36(/iK)2, /5sE = -3.0, QdE = 1.0(/iK)2, /3dE = 1-5, (7.6) 

QsB = 0.30(/iK)2, /5sB = -2.8, QdE = 1.0(/iK)2, /3dB = 1-5. (7.7) 

In both cases m = —0.6. Already from this rather naive argument (which can be applied, stricto sensu, 
only for sufficiently low multipoles, i.e. i < 100) interesting conclusions can be drawn. Consider, 
indeed, the ratio of the foregrounds: 

^^ ~ 0.5, z/ > 65 GHz. (7.8) 

7t(,ee} ^ / 

yi 

The frequency dependence of Eq. (17.80 (i.e. of the foregrounds) is much more shallow than the one 
expressed by Eq. (17.30 which would imply, in terms of the same ratio of Eq. (17.80 
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where u^ is an appropriate pivot frequency which can be chosen either to coincide with z/max or to 
be fixed experimentally on the basis of the specific features of the analysis. For instance, it could be 
useful to fix Vp to the value of the lowest (or the highest) frequency channel sensitive to polarization. 



^^The frequency slopes, directly known in limited regions of the spectrum, are often extrapolated over the whole 
physiccal range. 
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8 Magnetized CMB maps 

Satellite experiments like WMAP observe the cosmic microwave background in several frequency chan- 
nels. The observations are encoded in temperature and polarization maps. These have to be cleaned 
from foreground emissions, mainly coming from our own galaxy but also from extra-galactic sources. 
The resulting maps are used to deduce the angular power spectra C) introduced in Eq. (15.291) . 

Some strategies of parameter estimation are based on the comparison between the observed and 

(XY) 

the computed C) \ In a complementary perspective, the temperature maps already contain valuable 
information which can be used to distinguish between different cosmological paradigms such as the one 
suggested by the mACDM model. The local extrema in the temperature maps provide an interesting 
diagnostic that can be used to constrain the underlying cosmological scenario. Heeding observations, 
local extrema (i.e. hot and cold spots in the cosmic microwave background), have a large signal-to- 
noise ratio and are therefore easily detectable. On the theoretical side the statistics of local extrema 
assuming Gaussian fluctuations has been already studied in detail [72], [731 EH UHl [76] . 

To simulate CMB maps it is mandatory to have an efficient way of estimating power spectra in 
terms of a finite number of cosmological parameters. This step can be achieved via a Boltzmann solver 

(XY) 

such as COSMICS, CMBFAST and, among their descendants, MAGcmb. The calculated C^ can 
be fed back into the HEALPix package [771 EH] which provides a pixelization scheme for the data on 

(XY) 

the sphere. This is essential in order to pass from the angular power spectra Cg to a full sky map 
of the temperature fiuctuations and polarization as observed by different satellite experiments. 

Examples of sky maps will now be provided for different values of the pre-decopuling magnetic 
field and its spectral indejo. In particular, the attention will be focussed on three particular choices 
of parameters: 

• the best fit of the WMAP 5- year data alone (see Eq. (16. ip ): 

• the best fit of the WMAP 5-year data alone supplemented by a magnetic field with i?L = 5nG 
and spectral index ne = 1.5; 

• the best fit to the WMAP 5-year data alone supplemented by a magnetic field i?L = 20 nG and 
spectral index n^ = 2.49. 

The last choice of parameters is already ruled out by the observational bounds on the first acoustic 
peak, as discussed in Section [61 Nonetheless, it is instructive to consider also the last (extreme) 
model to emphasize a trend which is still there also for lower values of the magnetized background but 
which would be more difficult to visualize. The sky maps are simulated using the improved version of 
MAGcmb which has been described and illustrated in the previous sections. The magnetized angular 
power spectra have been calculated for a maximal multipole i^aax = 2500 and at a frequency 217 GHz. 
This frequency is the lowest Planck frequency with the highest sensitivity, at least according to the 
Planck bluebook [13] . Again, if lower frequency channels will be more sensitive than what claimed at 
the moment, this will be very interesting for our study, as the skilled reader can appreciate. 

For the simulations the value Ngide = 1024 has been adopted. This figure corresponds to a resolution 
parameter 10 used by the WMAP team. Furthermore, the first simulation was done using a Gaussian 



^^For the calculations reported in this section the use of the HEALPix package [77l [78] is warmly acknowledged. 
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beam with beam size of 5 arcmin, corresponding to the Planck channel at 217 GHz. In a second step the 
resulting maps were smoothed to a Gaussian beam size 6'fwhm = 1° i-e. of the order of the resolution 
of the five frequency channels of WMAP, which are between 0.22° in the highest frequency channel at 
94 GHz and 0.88° in the lowest frequncy channel at 23 GHz. In Fig. [12] examples of temperature maps 
at Planck resolution are shown. The effect of the presence of the magnetic field is an augmentation of 
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Figure 12: Temperature maps (in /xK) for the best fit WMAP 5 model without a magnetic field (left) 
and with a magnetic field B\^ = 20 nG and ne = 2.49 (right) simulated assuming a Gaussian beam 
size 6'fwhm = 5'. 

the hot spots as well as a general raise of their mean temperature as will be discussed below in more 
detail. The number of cold spots is increasing with increasing magnetic field strength and spectral 
index while their average temperature is decreasing. It can then be argued that the local extrema are 
becoming more pronounced as the magnetic field becomes more intense. A relatively strong magnetic 
field of i?L = 20 nG does change the appearance of the simulated temperature map. Conversely, a 
magnetic field of 5 nG does not have an effect simply visible in the temperature maps. The latter 
statement can be corroborated by looking at a 20° x 20° patch of the sky. Indeed, Fig. [13] shows, 
visually, that it is difficult to distinguish the imprint of a magnetic field strength i?L = 5 nG. By 
comparison, a 20 nG field can instead be distinguished (see Fig. [131 plot at the far right). The effects 
of a magnetic field can be also seen in the simulated polarization maps. In Fig. [14] the U polarization 
map is reported, just for purposes of illustration. The general visual result is that the U (and to a 
lesser extent the Q) maps are modified only for sufficiently intense values of the magnetic fields. What 
has been reported in Figs. [121 [IS and [H] represent clearly a feasibility proof of the strategy pursued 
in this paper. It also represents a novelty since, to the best of our knowledge, magnetized CMB maps 
have never been computed before within a well defined numerical approach, such as the one described 
here and initiated in [231 [211 [25]. 



In spite of the fact that CMB maps contain a lot of precious informations, they cannot be just 
analyzed with "visual" methods. We are going here to point out less anthropic ways of looking 
at magnetized CMB maps. A primordial magnetic field increases the number of extrema in the 
temperature map. Furthermore, it increases the mean temperature of the hot spots and lowers the 
mean temperature of the cold spots. To make this statement more prescise the HOTSPOT routine 
which is part of the HEALPix package has been used to find the extrema in the three cases at hand. 
This values of the local extrema and their localization can then be obtained. Suppose, for simplicity, 
we are only interested in the one-point statistics (hence the localization of the extrema is not taken into 
account). A similar approach has been taken, for instance, in [76j. The statistics of the temperature 
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Figure 13: 20° x 20° fields of tlie temperature maps (in /iK) in gnomonic projection for the best fit 
WMAP 5 model witliout a magnetic field (left), with a magnetic field Bi^ = 5 nG and riB = 1-5 (center) 
and with a magnetic field Sl = 20 nG and ub = 2.49 (right) simulated assuming a Gaussian beam 
size 6'fwhm = 5'. 
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Figure 14: U polarization maps (in fiK) for the best fit WMAP 5 model without a magnetic field (left) 
and with a magnetic field 5l = 20 nG and ne = 2.49 (right) simulated assuming a Gaussian beam 
size 6'fwhm = 5'. 

values (number, mean, variance, skewness and kurtosis) of the minima and maxima is illustrated in 
Table 1. For comparison, the simulations were done for A'^side = 1024 which results in A^pix = 12A'g^jjg 
given by Apix = 12582912. This means that about 4% of the total data points are local extrema. 
There are further quantities characterizing the statistics of the local extrema of the temperature map, 
like for example the two point correlation function. This has been used to investigate the detectability 
of weak lensing in the distribution of hot spots [75]. We leave this as well as other developments to 
forthcoming studies [82] . 

Smoothing the maps with a Gaussian beam of size 6'fwhm = 1° results in maps of lower resolution. 
This is equivalent to going from the resolution of the Planck experiment to WMAP. The resulting 
temperature maps are shown in Fig. [151 As it can be appreciated from Fig. [151 this resolution the 
effect of even a considerably strong magnetic field is hardly visible. 

Applying the HOTSPOT routine to the smoothed maps results in a considerably smaller number of 
local extrema. The trend of the magnetic field increasing the number of extrema is still observable. In 
this case, however, the increase is less pronounced, as can it be argued from Table 2. For comparison 
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Model 


type 


number 


mean(/iK) 


var(/iK)^ 


skew 


curt 


Bl = 


Max 


202174 


76.951 


5909.1 


-0.9012 


-0.0329 


5l = 5 nG, ub = 1.5 


Max 


203083 


77.778 


6000.0 


-0.9032 


-0.0340 


5l = 20 nG, ub = 2.49 


Max 


303089 


96.355 


7542.7 


-0.9728 


-0.0507 


Bl = 


Min 


201601 


-76.175 


5776.5 


0.9305 


-0.1108 


Sl = 5 nG, ub = 1.5 


Min 


202495 


-76.915 


5829.3 


0.9360 


-0.1093 


Bl = 20 nG, hb = 2.49 


Min 


302607 


-92.029 


6951.7 


1.0603 


-0.0949 



Table 1: Statistics of the temperature values (number, mean, variance, skewness and curtosis) of the 
local extrema (maxima/minima) of the simulated temperature maps at ^fwhm = 5' for the WMAP 
5-year bestfit model alone (i.e. Sl = 0) and for the WMAP 5-year bestfit supplemented by a magnetic 
field with strength i?L and spectral index ub- 



Smoothed mnP 5 alone beetf it, theta=l deg,20 x 20 (deg>"2 



Smoothed UMflP 5 + B=5rG, nB=1.5,th=l deg, 20 x 20 (deg) -2 



Smoothed mnP 5 + B=2QnG,nB=2.4S, th=l deg, 20 x 20 (deg>"2 






Figure 15: 20° x 20° fields of the temperature maps (in /iK) in gnomonic projection for the best fit 
WMAP 5 model without a magnetic field (left), with a magnetic field Sl = 5 nG and ^b = 1.5 (center) 
and with a magnetic field -Bl = 20 nG and ub = 2.49 (right) smoothed assuming a Gaussian beam 
size 6'fwhm = 1°- 



the statistics of the temperature values of the local extrema of the 5- year Internal Linear Combination 
(ILC) map provided by the WMAP collaboration has been included. This was obtained applying 
the HOTSPOT routine to the data file provided at the LAMBDA web site [79] . The ILC map 
is a weighted linear combination over the smoothed temperature maps obtained from each of the 5 
frequency channels. It has minimal Galactic foreground contribution and is assumed to give a reliable 
signal of the cosmic microwave background at angular scales greater than 10°. Its resolution parameter 
is 9, which is equivalent to A^side = 512 and it has a resolution of 1°. The number of local extrema in 
the simulated maps and the ILC map are of the same order. Comparing naively the mean temperature 
of the hot spots of the simulated maps with the ILC map seems to indicate that actually models with 
a non vanishing magnetic field are a better fit to the data. Likewise the models with a non vanishing 
magnetic field have a lower mean temperature of the cold spots which is in better accordance with 
ILC data then the WMAP 5-year best fit model without a magnetic field. To draw such conclusions 
one should use the WMAP 5-year temperature maps in the different frequency channels and perform 
a more thorough exploration of the parameter space of the models using Monte Carlo simulations |76j . 
This is beyond the scope of the current article [82]. The comparison with the ILC map was used for 
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Model 


type 


number 


mean(/iK) 


var(/iK)^ 


skew 


curt 


5l = 

5l = 5 nG, ub = 1.5 

5l = 20 nG, ne = 2.49 


Max 
Max 
Max 


7274 
7298 
7606 


63.005 

63.555 
66.198 


2715.5 
2731.2 
2868.1 


-1.0426 
-1.0484 
-0.1051 


-0.1742 
-0.1676 
-0.1843 


WMAP 5 ILC 


Max 


8066 


65.597 


2691.1 


-1.1370 


-1.1617 


5l = 

5l = 5 nG , ub = 1.5 

Sl = 20 nG, ne = 2.49 


Min 
Min 
Min 


7252 
7277 
7600 


-63.482 
-63.809 
-66.509 


2778.3 
2796.3 
2913.2 


0.9763 
0.9770 
0.9955 


-0.3041 
-0.3039 
-0.2846 


WMAP 5-year ILC 


Min 


8120 


-64.903 


2647.6 


1.1306 


-1.6289 



Table 2: Statistics of the temperature values (number, mean, variance, skewness and curtosis) of the 
local extrema (maxima/minima) of the simulated temperature maps smoothed with a gaussian beam 
with ^FWHM = 1° for the WMAP 5 bestfit model with and without a magnetic field with field strength 
-Bl and spectral index ub- For comparison the results of the Internal Linear Combination (ILC) map 
of WMAP 5 has been included. 



illustrative purposes showing that the local extrema of the simulated maps are good indicators of the 
one point statistics of the observed temperature map. 

Thus in conclusion, we see that primordial magnetic fields do have an effect on the structure of 
the temperature and polarization maps of the cosmic microwave background. These effects are clearly 
visible for rather large values of the primordial magnetic field at high resolution as expected to be 
provided by the Planck experiment. To quantify the effect of a magnetic field on the temperature 
map the statistics of the local extrema has been used. On the basis of the reported results, the Planck 
experiment provides an interesting perspective of constraining magnetized models using temperature 
maps. 

9 Concluding considerations 

In the present paper a viable strategy for the complete calculation of the magnetized polarization 
observables has been presented. To achieve this goal it was mandatory to employ, within the different 
branches of the spectrum of plasma excitations, the appropriate description. In previous studies 
this aspect was never discussed in detail and never applied. The reported results are a necessary 
physical step for any consistent calculation of the magnetized polarization observables either at the 
semi-analytical or at the fully numerical level. 

At a more operational level, the accurate calculation of the TE, EE and BB angular power spectra 
beyond the MHD approximation, will be of upmost importance for our program of parameter estima- 
tion which is already under way. In front of difficult problems, such as the one we are scrutinizing, 
there could be the tendency of substituting words for actions. Even if this tendency could be poten- 
tially rewarding in the short run, it will be the bane for future comparisons with observational data. 
Our primary objective, in this respect, is very modest: we want to develop a sound approach where 
all the relevant physical effects (i.e. dispersive and not dispersive) are realistically modeled. Absent 
this necessary step, any parameter estimation would be forlorn and theoretically biased. 
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The second line of developments of the present considerations imply a strategy for deducing the 
typical signatures of a Faraday-induced B-mode. The reported results suggests that an unambiguous 
way of pinning down the Faraday-induced B-mode is to look at the scaling properties of the measured 
BB angular power spectrum as a function of the frequency and to compare both with computed BB 
power spectrum as well as with the other power spectra of the cross-correlations (i.e. EE and TT 
angular power spectra). 

Realistic simulations of maps including the effects of large-scale magnetic fields have been derived 
and presented. For experiments with high (nominal) resolution (such as the one of the forthcoming 
Planck explorer mission) the effect of a primordial magnetic field on the temperature and polarization 
maps is distinguishable. Lowering a bit the resolution changes the structure of the temperature maps 
but, according to our preliminary results, does not totally jeopardize the possibility of reading-off the 
effects of (sufficiently strong) pre-decoupling magnetic fields. Polarization maps have been also studied 
for different magnetic field parameters. For purposes of presentation, only few examples have been 
illustrated. As a general rule of thumb, a magnetic field increases the number of hot and cold spots in 
the simulated temperature maps. This has been quantified in the determination of local extrema and 
their one-point statistics. There is a potential correlation between the number of local extrema and 
their corresponding temperatures. The differences between the mean temperatures for the hot and 
cold spots, respectively, for different values of the magnetic field strength and its spectral index are 
larger at a higher resolution of the experiment, or in other words for a smaller beam size. Our results 
open interesting perspectives for the Planck explorer mission: the mACDM scenario (as well as any 
of its non-minimal extensions) can be analyzed by studying the one- and two-point statistics of the 
temperature maps 

The original spirit of our endeavors has been, also in the past, to bring the treatment of magnetized 
CMB anisotropics to the same standards typical of the conventional case. At the moment, a number 
of preparatory tasks for the forthcoming Planck experiment has been completed. The fulfilled tasks 
permit, at the moment, explicit and accurate calculations of the effects of pre-decoupling magnetic 
fields on the various CMB observables within a faithful dynamical framework. 

There are other intriguing physical problems which are linked to the ones we discussed in this 
paper. An example along this direction is already apparent from the scaling properties (with frequency) 
of the temperature and polarization autocorrelations which has been suggested, in this paper, as a 
model-independent diagnostic of magnetized birefringence. It is here important to stress that the 
experimenters will have to apply a mask selecting those regions with low contamination of Faraday 
rotation due not to the pre-decoupling magnetic field but rather to the galactic magnetic field. It 
seems that here there will be, in the months to come, a rather interesting interplay between the latter 
problem and the recent findings of cosmic ray experiments lj as far as the properties of the large-scale 
magnetic fields in the local Universe are concerned. This problem, for the reasons mentioned above, 
is also rather interesting for us [82] . 



■^^The latest analyses of the AUGER experiment demonstrated a correlation between the arrival directions of cosmic 
rays with energy above 6 x 10^® eV and the positions of active galactic nuclei within 75 Mpc |£0i|. At smaller energies it 
has been convincingly demonstrated [5T] that overdensities on windows of 5 deg radius (and for energies 10^^'^eV < E < 
10^^'^eV) are compatible with an isotropic distribution. Thus, in the highest energy domain (i.e. energies larger than 60 
EeV) , cosmic rays are not appreciably deflected: within a cocoon of 70 Mpc the intensity of the (uniform) component 
of the putative magnetic field should be smaller than the nG. On a theoretical ground, the existence of much larger 
magnetic fields (i.e. 0(/iG)) cannot be justified already if the correlation scale is of the order of 20 Mpc. 
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A Vlasov-Landau equation in curved space-time 

The Vlasov-Landau equation for charged species is rarely written in the presence of relativistic fluctu- 
ations of the geometry. In this Appendix the Vlasov-Landau approach [83l [84] (see also [85j) will be 
derived in curved space-time l1| for a generic species which we will take to be positively charged. The 
extension to negatively charged species is trivial. 

In a curved background, the best pivot variables for the Vlasov-Landau approach are the comoving 
three-momenta (very much as it happens in the case of the conventional Boltzmann equation which 
is implemented in standard Boltzmann solvers). In the present case the main difference is that the 
geodesic of a charged species is affected by the electromagnetic fields. Consider, to begin with, the 
conjugate momenta P" and their mass-shell condition for a generic (massive) species: 

(iT 

P° = mw" = m-— , gapP^P^ = m^, (A.l) 

where m" is the four-velocity of the generic species and where A denotes, throughout this Appendix, 
the affine parameter. The mass-shell condition of Eq. (lA.ll) implies, as it has to, gapu'^u^ = 1- 
From the second relation of Eq. ( lA.ll) the three-momentum p* can be introduced in terms of the 
conjugate momenta as gijP^P^ = —Sijp^p' . The comoving three- momentum is then, by definition, 
g* = ap*. According to Eq. ( JA.ll) the components of the conjugate momenta and the components of 
the comoving three-momenta are related as: 

P° = \jm^a^ + q^, Pq = Jm^a^ + q^, (A.2) 



a2 



q' = a' 



■ K 



pi, (A,3) 



where hij is the metric fluctuation defined in Eq. (12. 2p . It is useful to relate the comoving three- 
momentum related to the comoving three- velocity v which has been used in the bulk of the paper and 
which arises naturally in the kinetic treatment. The relation between v and u is given by: 

u = -7- = M f , « = IT = ^-^) 7(^) = . (A.4) 



dX dX a a/1 — 



V 



2 



Since, by definition, f* = P''/P^, we shall also have 



q _^ _^ 

q = ma'y{v)v. (A. 5) 



7' 



Having introduced all the relevant variables we can write the geodesic for charged particles, namel}o: 

^ + r>"«^ = -F'^"w„, (A.6) 

aA '^ m 



which becomes 



dv V e , ^ ^, . . X 

:r + ^- = — ?^ + ^ X p , A.7 

dr 7^ may^ 



dr mc 

where E = a?£ and B = a?B. Concerning Eqs. (jA.TD and (JA.SP two comments are in order: 



^ + Hu\l + v^) = ^E-v, (A.8) 



^^Even dedicated discussions of kinetic theory in curved space-times [S^, do not address the problem of the Vlasov- 
Landau description in curved backgrounds when curvature inhomogeneities are simultaneously present. This is one of 
the aims of the present Appendix. 

■^^We will consider, for sake of concreteness, the case of a particle with charge -^-e and mass m. 
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even if we used rescaled electromagnetic fields, the evolution equations are not the ones we would 
have in flat space-time; 



Eq. ( 1A.8|) is implied by Eq. ( 1A.7I) once we recall that, by definition, m° = 7(f) /a. 



The rationale for the first comment is that, of course, the particles are massive, in our case and the 
very presence of the mass breaks Weyl invariance. Notice, however, that by using the definition of the 
comoving three-momentum g* in terms of the comoving three-velocity v (i.e. Eq. flA.5|) ). Eq. flA.7l) 
can be written as 

^ = e(E + vxB). (A.9) 

dr 

Now, in the ultrarelativistic limit v = q/\q\ (and the evolution equations would have the same flat- 
space-time form). In the non-relativistic limit v = q/{ma) (and conformal invariance is now broken). 
Clearly, in the problem addressed in the present paper the electrons and the ions are always non- 
relativistic throughout all the stages of the calculation. Recalling that, to first order in the synchronous 
fluctuations of the geometry, 

r?, = ns,, - \{K^ + 2'Hh.,), r|o = ml - h{, (a.io) 



Eq. (jA.OP gets modified as 



■mav^'-f{v) 






(A.ii; 



The Vlasov-Landau equation can then be written as: 



dl idl_ 
dr dx' 



-j-q^ + e{E' + e'^^'vkBj) 



^ = CcoU. (A. 12) 



By taking the integral of the right hand and of the left hand sides of Eq. (1A.12P we do obtain the 
canonical form of the continuity equation 

On -> h' 

— + V-{nv)--n = 0, (A.13) 



where we assumed that the massive species are non-relativistic. The result of Eq. flA.13|) can be easily 



obtained using integration by parts exactly as in the flat-space analog [T]. Notice, indeed, that the 
integral of the collision term over the velocity vanishes and that the derivative of the Lorentz force upon 
the velocity also vanishes. Equation ( JA.131) can be obtained directly from the covariant conservation 
of the current, i.e. 

^t.r = 0, r = efiu^ (A. 14) 

by recalling that V^j^ = d^j^ + r^^^j" and that n = a^n. Neglecting hij Eq. flA.12p can be written, 
in the non-relativistic limit, as: 

^ + v% + —{E + vxB)% = Coll. (A.15) 

OT OX ma ov 

Equation flA.lSp describes, for instance, the evolution of the electrons and ions prior to equality and 
throughout decoupling. In this case the equilibrium distribution will be Maxwellian. At the same 
time Eq. flA.121) has been deduced without specifying the equilibrium distribution and it can then be 
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applied in more general terms. Consider, indeed, the situation of massless neutrinos. From Eq. (]A.12p 
setting E = and B = we get 



df_ (h^dl_ (fq^df_ ^ 

dr dr dx^ '^ 2q dq 



(A.16) 



For a massless particle gapP'^P^ = and x^ = rn^ and g* = n^q. Thus, Eq. (JA.161) can be written as: 

r)f f)f 1 rif 

(A.17) 






By perturbing Eq. (1A.17P to first order in the metric inhomogeneities (around a Fermi-Dirac distribu- 
tion f^^\q)) we obtain: 



dfW 
dr 



+ ikfif 



(1) 



r^^^ 



ih' + 60 - r 



dlnq 



0, 



(A.18) 



where fi = k-h and the fluctuation of / has been defined as f{q, k, fi, r) = f^^\q)[l + f^^\k, jj,, t)]. By 
now defining the reduced phase space distribution for the neutrinos, i.e. 



Tu{k,^i,T) 



Iq'dqf('\q)f('\k,fi,T) 
jq^f^^){q)dq 



(A.19) 



we can easily obtain the result reported in Eq. (I3.25p . The quantity J?v(/c, /i,r) is the one usually 
expanded in Rayleigh series. 



^,(A;,/i,r) = 5:(-0^^,,(A;,r)P,(/i), 



(A.20) 



where Tvt, are the multipole moments. In the present discussion, the initial conditions of the neutrino 
hierarchy are set by solving for the lowest multipoles, i.e. monopole, dipole and quadrupole 



B Magnetized adiabatic mode 



As discussed in the bulk of the paper, the initial conditions for the full Boltzmann hierarchy are given 
well before equality, i.e. in a regime where the ions, the electrons and the photons are all coupled 
together for different physical reasons. Given the largeness of the Coulomb and Thompson rates in 
Hubble units, it is practical to set initial conditions for the MHD fluid (i.e. the baryon fluid in the 
conventional terminology). The structure of the magnetized adiabatic mode has been derived, in the 
synchronous gauge, in Ref . [H] . Previous analyses in different gauges can be found in [T^l [201 EI] 
and also in [ST] (also in the synchronous gauge). The magnetized adiabatic mode is here reported for 
completeness: 



e(fc,r) 

h{k,T) 



-2C{k) + 



4i?, + 5 



6(4/?, + 15) 



C{k) + 



R^iAanik) - R^n^ik)) 



-C(k)k'T' - — 
^ ' 36 

R^{15-20R, 



^Rl - lAR,^ 



75 



6(4i?, + 15) 
-C{k) 



k^T\ 



(B.l) 



10(4i?, + 15)(2i?, + 25) 



(2/?, + 25)(4i?, + 15) 



4^4 



A;V' 



(B.2) 
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S,y{k,T 

6c{k,T 

Sb{k,T 

9-yb{k,T 

e,{k,T 

ec{k,T 

a^{k,T 



-IL,nBik) 
-R^nBik) 



C{k)-an{k) + ^nn{k) 



2_2 



k'r 



C{k) + 



R-, 



^Ru 



Aa^ik) - R^n^ik) 



2_2 



Fr 



--R^n^^k)- 2 



C(A;)^2^2 



--R^n^ik) - 



R 



C{k)-aB{k) + ^nB{k) 



-^Vljiik) -cxb 



36 



2C{k) + 



fcV, 



4^3 



rr 



IjaB(fc) - ^1^b(A:) 



fcV- 



36 



2(4i?, + 23) 
AR^ + 15 



C(fc) 



+ 



2R,{AR, + 15) 
0, 

It,, 



(AaBik) - R^n^ik)) 



4^3 



A;V 



4C(A;) 



+ 



/?^(4aB(A;)-i?,f}B) 



3(4i?, + 15) 2R,{4R, + 15) 



2_2 



k'r 



(B.3) 
(B.4) 

(B.5) 
(B.6) 
(B.7) 

(B.8) 

(B.9) 

(B.IO) 



Note that, as thoroughly discussed in [241 125] the curvature fluctuations on comoving orthogonal 
hypersurfaces can be simply expressed in terms of ^ as 



7^ = ^ 



m' 



The power spectra of curvature inhomogeneities will be assigned as: 

27r2 



(7^(A;)7^(p-)) 



A;3 



■Vn{k)6'^''\k-p). 



(B.ll) 



(B.12) 



i.e. in full agreement with the way the magnetic power spectra have been assigned in Eq. (13.291) . The 
two-point function of Eq. fl3.29p is completely specified only if the magnetic power spectrum is given 
in explicit terms. We will choose both for the power spectrum of curvature perturbations and for the 
magnetic power spectrum the following power-law parametrization: 



h \ ris — 1 

Vn{k)=An[Tr 



Vn{k)=AB{- 



k \"B- 



(B.13) 



where fcp and /cl are the two pivot scales, rig and ub are the two spectral indices. In the conventional 
ACDM scenario the amplitude of curvature perturbations is customarily given in terms of A-r which 
is,by definition, the amplitude of the power spectrum evaluated at the pivot scale kp. In the case of 
the magnetic fields we will assign not the amplitude of the magnetic power spectrum but rather the 
regularized field. This technical aspect will be recalled in the following Appendix since it is required 
in the estimate of the angular power spectrum of Faraday rotation. 



C Power spectra of Faraday rotation 

If the magnetic field is not homogeneous, the Faraday Rotation rate will be, itself, not a homogenous 
quantity [16]. Here the power spectrum will be computed. This calculation was also discussed before 
in the literature (see, in particular, [39] and also [IQ]). The notations employed in the present paper 



44 



are slightly different due to the fact that we are presenting here a full calculation of the magnetized 
polarization observables. The Faraday rotation rate can be usefully expanded in ordinary (i.e. scalar) 
spherical harmonics: 

F{n) = J2fen.Yem{n). (C.l) 

im 

Since the geometry is isotropic the ensemble average of the fim must not depend upon m but only 
upon i and, consequently, on symmetry ground we can expect that 

{flmfi'm') = Cl Sei'Smm'- (C.2) 

By using the addition theorem of ordinary spherical harmonics the two-point function of the Faraday 
rotation rate can be written as 

(F(n)F(m)) = ^ ^^^C^ P^ih ■ m). (C.3) 

Equations (lC.2p and (1C.3P are only meaningful if Cf is computed in terms of the magnetic power 
spectrum. To achieve this goal it is practical to expand the magnetic field not in scalar spherical 
harmonics but rather in vector spherical harmonics [891 EH]- In the latter framework it is possible 
to find the vector analog of the well known Rayleigh expansion. As thoroughly discussed in [88j this 
technique is rather well established both in the study of electromagnetic processes [SE] as well as in 
nuclear physics \UU\. The vector analog of the Rayleigh expansion can be written, in the present case, 
as: 

me-'''^^' =EY.9itik,^^)Y}:iih), (C.4) 

where, as usual jj, = k-n. In Eq. (IC.4p a is the polarization and F/^ are the vector spherical harmonics 
or, for short, vector harmonics in what follows. They can be written as 



Yt!(n)^^^M, y-l«)(fl) = -iiii^Z^ffi^. ?/.;■'= ny,„(fi). (C,5) 



+ 1) l/C(C+l) 

where Yem{n) are the usual spherical harmonics. Sometimes 1^ J(n) and F/J(ri) are called, respectively, 
the electric and the magnetic multipoles. We shall avoid this terminology which could be potentially 
confusing. 

Equation flC.4p holds for a generic vector, we do know that the magnetic field is transverse and, 
therefore, the expansion coefficients are more than needed. Furthermore, it should be borne in mind 
that, in the Faraday rotation rate F(n) the magnetic field appears to be projected over the direction 
n. The two previous observations immediately demand that: 

n ■ F£(n) = n ■ Y,^'2in) = 0, k- B{k) = 0. (C.6) 

As a consequence of Eq. (]C.6p . the multipolar expansion of Eq. (1C.4P simplifies: 

h- B{k)e-''^^^ = Y.gUk.^^)n■Yi;,'\n), (C.7) 

where gim{k,fi) in Eq. (IC.Tp is the only non-vanishing coefficient of the multipolar expansion and it 
is given by [89] : 

9im{k,fi) = ^""27+1 \ ^i{i+l)[je+i{kTo) +je^,{kro)]B{k) ■ Y}'2*{k) 

- [{£ + l)j,-,i(A;ro) - ije-i{kr,)]B{k) ■ Yt'^*{k)}. (C.8) 
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Once more, the magnetic field is transverse (see, e.g., Eq. (1C.6P ). Thus, since 1^ ~ (k) is proportional 
to k, the last term in Eq. (IC.Sp is identically zero. Recalling the well known recurrence relations of 
(spherical) Bessel functions [91], [92] , 

2i + 1 

je+i{kTo) +je-i{kro) = — ji{kTo), (C.9) 

fcro 

Eq. (1C.7I) becomes: 



n ■ 5(fc)e-^^^- = 4vr ^ ^i{i + lY-^B{k) ■ F/^^ * (k). (CIO) 

From the result of Eq. (JC.IOI) the f£m of Eqs. (JC.ip and (10.20 can be swiftly determined and they are: 

/.^ = ^^§^^^fii^^)Jd'k'^m-^ (c.ii) 

Using Eq. (13.291) together with the orthonormality condition of vector harmonics, i.e. 

J dnf^Y}^{k)Yl^>{k) = 6u'Smm', (0.12) 

the angular power spectrum of Faraday rotation can be simply expressed as: 

CP = AnAH{l + 1) / ^Vsikf^^j^ (0.13) 

Since, by definition. 



I TT / k \ "■B-1 

MkTo) = J^J^+Mkro), VB{k) = ^b(^J , (0.14) 



Eq. (10.13P can also be written as: 






CY> = 2'k'(.{(. + 1)A'Ab[^] J(nB,£), (0.15) 



-pl 5-nB \-p\ P I "B 3 

/•OO 1 \ 9 I \ ~'~ 9 9 

^^K,^) = / xl-~'Jl,/,{x,) = -^^ f^ f, (0.16) 

JO /VTT -p/ 6-nB \r I 7 I n ns 



2V^y(^\v{i + 1- 



where we used that xq = kr^ and that fco = Tq^. Note that the integral of Eq. (10.161) converges for 
— 1 < riB < 5 so there is no need of an ultraviolet cut-off o To make the expressions even more 
explicit we will adopt the following strategy: 

• first we will trade A^ for the regularized magnetic field intesity -Bl; 

• second we will refer the frequency of the channel u to the frequency of the maximum of the 0MB 
emission. 



■^•^This is also what happens, for instance, in the case of the exphcit calculation of the Sachs- Wolfe plateau where a 
cut-off is not usually imposed because the integral is convergent in the physical range of spectral indices. 
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Complying with the first step we have that: 



(2vr)"B-^^2 



r(^ 



L5 



riB > 1, 



A, 



l-UBfkoY-^^ . 



h 



L5 



TiB < 1. 



(C.17) 



(C.18) 



Equation (]C.17p adopting a Gaussian window function, while Eq. (JC.ISP is derived by using a step 
function regularization cutting, effectively, all the modes with k < ko [2^. 

The maximum of the CMB emission can be easily derived by maximizing the energy spectrum of 
the CMB, i.e. 






15 x^ 

r' '7O" 



fi^o = 2.47 X 10~^ 



where x = k/T^Q. By taking the derivative of Q^{x) we can find the maximum: 



dx TT 

From Eq. (lC.20p we do have that 



15 a;3[4e^(x-4)] _ 

~~^^'"tO ]— TT — 0; 



X-n 



3.920 






27r 



'70 



222.617 GHz. 



Consequently, the final formulas for (7) ^ can be expressed as: 

Cf^ = CZP{nB), nB>l, 
CP=Czf\nB), nB<l, 



^(F) 



(F), 



where C, Z) [n^ and Zn [n^ turn out to be: 



c^5£!e.-n 



lOe' 



2 max 



BL 



V„ 



30.03 ^BL 



V 



v„ 



^ 



_ Bl 

BL — o - ' 

87rp^ 



;r(F). ._(ko\---HU + 



i)(2.rB-r(^)r(£ + ir-|) 



zf^K) 



1 -riB 



rf^ r 6^rU + £-f 



+ i)r(^ r £ + ^-| 



r(^ rU + 



2 



(C.19) 



(C.20) 



(C.21) 



(C.22) 
(C.23) 



(C.24) 
(C.25) 

(C.26) 



where: 



(BB) _ T^j2 Y- ^r2 Vie I) I) ^2( 2^1 + 1)(2^2 + l) ^(EE)^(F)r^ffl n2 

47r(2£ + 1) ^^ ^^ ^ ^^°^^°^ ' 



cr^ = iVi^E^i.^(^,^i,^2 



^1,^2 



/C(£, £1, £2) = --[iv' + L^ + L^ - 2L1L2 - 2LiL + 2Li - 2L2 - 2L], 

(^2 + 1)^2, 



L = (£+l)£, Li = (4 + l)£i, L2 


Nn 


2(^-2)! ., 


2(^2 -2)! 


J-^e. — \ 


(£ + 2)!' ^ '^ \ 


(4 + 2)! 
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(C.27) 

(C.28) 
(C.29) 

(C.30) 



In Eq. (1028]) Cfo^^o is given by: 



Wl0f2 



'-'£10^20 



+ ii + 



2n + l, 



(n-£i)!(n-£2)!(n-£)!\ 



(2n - 2£i)!(2n - 2£2)!(2n - 2£)! 
(2n + l)! ^ 



+ £1 + 



2n, 



where n is a positive integer. This form of the relevant Clebsch- Gordon coefficient is given by [89] (see, 
in particular, page 251) and it has been also used in [IQ] where Eq. f lC.27p has been firstly derived. 
The Clebsch- Gordon coefficient of the previous equation then vanishes unless ji'i — £2! < i < h + (^2 
(triangle inequahty) and unless ii+i2+^ is an even integer. In the two degenerate cases (i.e. i = £1+^2 
and ^ = ^1 — ^2) the expressions become, respectively: 



/^^l+£2 
'-'£10^20 



pil-i2 

'-'£i0£2 



+ 



(2£i)!(2£2)! 

iiU2\ \{2i, + 2i2)\ 

-1 



(2£i)!(2£i- 2^2 + 1 



4!(4-^2)!\ (2£i + l)! 



It should be borne in mind that, for some applications, it useful to deal directly with the amplitude 
of the magnetic power spectrum A-q- Following the conventions expressed by Eq. 03.29p . A-q has 
dimensions of an energy density. In terms of A-q we shall then have 



C 



Co 



Coi{i + l) 



r(^ r £ + ^-| 



n^iri, 



2 



A;o\"B-i 



2.863 X 10 



-6 



^-' 



nG 



u 



(C.31) 



In terms of the parametrization introduced in Eqs. (I3.29P and (IB. 131) it is clear that As has exactly 
the dimensions of a magnetic energy density. 
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